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ANALYIS OF PROPERTIES AND SYNTHESIS OF AUTOMATIC 
CONTROL SYSTEMS WITH LAG 


Fan Chun Wui 


(Moscow) 


A method of synthesizing compensating devices of servo sys- 
tems [1-3] is extended to handle systems with lags. Curves are de- 
veloped for synthesizing systems with lags and for analyzing their 
properties, Circle diagrams are constructed for the determination 
of the distortion function of the real frequency characteristics for 
different values of lag. 


In most works concerned with automatic control systems containing lags, attention {is principally directed 
towards the investigation of stability [4, 5]. 


In the present paper a method {s presented for estimating the influence of lag on the properties of an 
automatic control system, a method which makes possible the synthesis of the proper compensating device by 
using the frequency characteristics of the system which would obtain were the lags not present. 


1. Initial Relationships 





It will be assumed that the stabilizing feedback path does not shunt any lag elements, and that the system 
under consideration is stable for sufficiently large gain [6]. 


We denote by KW(s) the transfer function of the open-loop system, excluding the lag element. Then, the 
transfer function of the closed system, having but one feedback path, will have the form; 





_ KW(s)e—** 
pee rg 1+ KW (s)e—**’ (1) 


where r is the lag time. 


We introduce the following notation: 


KW (jw) = U (w) + jV @), O.(jo) = P;(@) + 70s (@). 


Substituting s = jw in Expression (1), we obtain the system's frequency characteristic; 





. U “pe 
Selle) = Be tare othe ; 


We shall consider the real part of , (jw): 
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a U?+V2+ U cos ta +-V sin ta 
+ 14+ U2 4 V2 + 2(U cos ta + V sin Tw) * (3) 





We introduce more notation: 








4 
F = U (cos tw — 1) + V sin tw a 
(i+ U)+V? : 
G= ars U (cos to — 1) + V sin ta (5) 
~ 14 2F ~ (14+ U0)* + V2 + 20 (cos tw —1) + 2V sin tw ’ 

— Re KM Lie) _ 

P= Re iT ewge 
(6) 

After a simple transformation, Expression (3) may be written in the form 
i 

Pe =a ar = P—G(2P —1). (7) 


As will become clear in the sequel, G (w, r ) of (7) can be considered as the distortion function. 


Expression (7) gives the relationship between the real frequency characteristics of the system with, and 
without, lag. It is known that the real frequency characteristic completely determines the behavior of the 
system with a single step-wise action, and has also some bearing on the behavior of systems with other modes 
of action. The problem thus reduces to the study of the influence of lagging links on the nature of the varia- 
tion of the real frequency characteristic. We shall assume that the form of the function P( w) has been chosen 
to meet the qualitative conditions affecting the synthesis of the system when lag is not taken into account. In 
this case, the term G(2P-1) completely depends on the function G. 


The problem thus consists of the following: 


1) to find the relationship between the distortion function G(w, 7) and the parameter r , i.e., to find the 
quantitative change in the real frequency characteristic when a lag of time r is present; 


2) for a given value of r , to carry out the synthesis of a compensating device which will allow the given 
qualitative demands on the system to be met. 


Consider the function G(w, 1). After a slight transformation, we obtain from Expression (5) 








(2G—1)costw+i 7? (2G—1)sintw }? (2G — 1)(cos tw — 1) 8 
ee te eC 


For a given value, G = const; by changing the value of rw we obtain the equations for a family of circles 
in the KWjw) plane. Holding rw = const and varying G, we obtain another family of circles (Fig. 1-5). 


2. Circle Diagrams for the Distortion Function Gw, r). 





We wish to dwell now on a question which has great practical interest. Let us assume that G = Ga and 
TW = gg. Substituting these values in (8) we obtain a circle, A, on the KWjw) plane (Fig. 6). On the same 
plane we plot the frequency characteristic, KW(jw). We assume that the circle A and the frequency charac- 
teristic KW( jw) intersect in the points a and b, corresponding to the frequencies wa and wp. On the basis of 
the circle diagrams, Fig. 1-5, we have that, if the lag time of the system r=1T, = ¢, /wg then, in the fre- 
quency interval 0 < w < wg the very important inequality holds; 


|G|<|G,. (9) 
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Fig. 1. Circle diagram for G = —0.1 and various Fig. 2. The same as Fig. 1, but magnified four times. 
values of rw. 











Fig. 3. Circle diagrams for G = 0,1. Fig. 4. Circle diagram for r w = 30° and various values 
of G. 


This inequality does not hold at point b, For r = gg / wh in the frequency interval wa < w < uw, 1G, may 


be larger than |Ga, since all the remaining circles for |G, > |G, andr w < gz lie inside circle A. We shall 
call the point a the first point of intersection, and point b the second, 
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Fig. 5. Circle diagrams for rw = 10°. 








We consider the circle diagram for G = — 0.1 
(Fig. 7). On the KW( jw) plane we construct the charac- 
teristic, KW(jw), which is intersected at point C, cor- 
responding to the cut-off frequency we, by one circle 
of the proper radius. Through point C there passes 
just one circle of the family, G = const. We assume 
that point C is the first point of intersection, From 
this follows the obvious result; if the lag time in the 
system TS Tco= ¥ [We (in the given case, the angle 
Yo equals 6°) then, in this system with frequency w< we, 
the value of |G|can not exceed 0.1. Thus, the lower 
i$ We, the higher is the maximum value ofr necessary 
to satisy Condition (9). The function G determines the 
influence of the lag. Hence, by use of the diagram for 
the given system, one can always find an admissible 
T max Which will guarantee that the actual distortion 
of the real frequency in the interval 0=w=w, does not 
exceed a previously given maximum value of |G |. 
Such a limitation on lc| ‘can be made in the widest in- 
tervals, for example in 0s w=<w,. Then, instead ofthe 
circle P, = 0.5, we choose the circle P,, = 0, and point 
II determines some value of t max- 


We consider now the obverse problem. If KW(jw) 
and r are given, then one may find the maximum value 
of the distortion coefficient, |G | , in, for example, the 
interval 0<w=w,. For the solution of this problem, 
one may make use of Table 1. 














jV 
nw( jus) Table 1 is entered with specific values of ¥, = 
=" +0 (wc) and g= wr. The quantity y, can al- 
g v ways be found from a known equation if KW( jw) is 
known: 
@b 
a 
Y= t+ 4(e) = (10) 
aio, 
ee ee ae 
@ - m—V¥>— >j tan? Pie + Si are tg Tice. 
i=] =] 
Fig. 6, 
TABLE 1 
yedeg. | 10 | 2] 9] 40 | 50 | 6 | 7 w |» we | 110 120 | G, 
0.5)1 {|1.5] 2 2.5] 3 3.5) 4 4.7| 5.2] 6.2] 7.4] —0.05 
1 2 i3 4 5 6 7.2| 8,7 9.8 | 14.2 | 13 15 —0.10 
1.51}2.814 5.3] 6.8] 8.4] 9.8] 11,8 | 13.3 | 16 18 20.5} —0.15 
?, deg. 2 {4 115.6] 7 9 /41 13 15,5 | 17.5 | 19.5 | 23 26 —0.20 
2.515 7,9}10 12.5/17.5| 20 22 25.5 | 30 35 49 —0.30 
$.2/5.5/8.5/11.2/)13.5/18.2| 21.5] 24 27.5 | 32 37 41 —0.40 
4 |6 /|9.5)12.5/16 19 |23 29 32.5 | 36.5 | 40 44 —0.50 
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The quantity wer is known if the system lag time, 7, is given, From these parameters, it is possible 
to find G, from Table 1. If point C is the first point of intersection, then the value of G, is indeed |G ahd 
in the given interval. 


We now consider the case when point C of Fig. 7 is not the first point of intersection (Fig. 8). The fre- 
quency characteristic first intersects the circle at point C, and thereafter at point Cg. In this case, the cut-off 
frequency corresponds to point Cg. In order to satisfy the inequality |G |< |Ga| for 0 =w = We, it is nec- 
essary that r be limited by the inequality r= 9 /we instead of r= ¢¢q Awe, where yogis the index of the 
circle tangent to KW(j w) at some point Cy (in the given, gcg= 4°, cg = 6°). 


In first and second order static, and astatic, systems, G can not be positive for low frequencies, but for 


astatic systems of the third order, the low frequency portion of the KW(jw) curve is intersected by circles cor- 
responding to positive values of G (Fig. 3). 


3. Method of Analyzing, by Means of the Distortion Function G, the Influence 
of Lag r on the System's Performance Quality 








It is well known that, by the principle of superposition in linear systems, the influence of a lag element 
on the transient response of such a system with step-wise response may be determined immediately from the 
second component of Equation (7). Denoting this term by ¢€ (t), we have 


Z, (t) = x(t) + e(t), 


where x, (t) is the transient response taking lag into account, x(t) is the transfer response without lag being 
taken into account and 


e(t) = — =| G(2P—1) 2 da (11) 
0 


@ 





We mention some properties of the distortion function G. 


It is clear from Formula (5) that for sufficiently small values of r w, the function G may be written in 
the form 
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Vt 
C= ayo TT = Om, (12) 





where 


he KW (jo) 
c= im 7 KW (jo) * 


In the majority of cases for ordinary systems, {Q| » for low frequencies, increases linearly, i.e., 























Qs < notin w. This means that, for low frequencies, G varies according to the law 
d 
G = we zw? (13) 
( cf, curve 2 on Fig. 9). 
-¢ oe ad 
op heaton Anax 
] } 
| 
| 
1 | Qs 
, 
le . 
a “ 
Fig. 9. Fig. 10. 


If it is known that, for w < Wp: |G | is always less than a certain value, |Gp | then, for w<wp, the 
form of the function G is approximately as shown in Fig. 9 (curve 2). 


In order to simplify the subsequent discussion without thereby introducing significant error, we shall 
assume, on the basis of the representation given above, that in the interval 0=w= wp G varies linearly, with 
slope given by tan a = Gp / wp, ie., 


G 
one ( for © S ®)). (14) 


We now find an approximate value for the second member of the transient function ¢ (t). It will be 
assumed that the real frequency characteristic, P (w), is of the type shown in Fig. 10, the parameters of which, 
as determined by the conditions on its properties, may be the following [1): 


2 @ 
hams <8, hae > 0.5. X, = (15) 


@, 

®p 
 & > 0.4. 

On the basis of Equation (11), and taking Condition (14) into account, one may write 
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On 





2¢6 in wt 
¢(t) =— = | 2o(2P—1) dw— 
9 2 ‘ (16) 
— ) G(2p—1) de, 
“p 


Initially, the second member is negligible. As in the analysis of an ordinary system (Fig. 11), its value 
may be estimated from the following inequality [1}: 

















2 C i 
|e, (t)| = =\ — G(2P — 1) “2% deo|< 0.3|— G (2P — 1)min| 
“P 
or “ 
|ta(t)| =| = (P +1—G(2P — 1) dal < 

“Pp 

<0.3\{P + {| — G(2P— 1)}}min. (17) 

Moreover, it is known that the entire real fre- 

quency characteristic can be represented as a collec- 
Di tion of a certain number of trapezoidal frequency 
«) = 4 characteristics, and (16) ean be written in the form 
\ 
; 2 —1 
Y w — ) = a > ( sin wt dw +- oes a , (18) 
0 S77 P pT? ”P 
Fig. 11. 


where P; is an elementary trapezoid of the collection. 


b 
eft a E(t 
. t Pp =f? fe 
ma 
a4 a¢ 
02 wot @ 





O x 2n inda Sx in mt bt on 


r -# 


Ox Dy 9% 4% Sa in x bx 





2 
i 6 ae 
04 
az m 
Ox in in 4n Sn ine bt ot O «© oy in 42 5x bt Tx bn Ox 
Fig. 12, 
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For the typical real frequency characteristic (Fig. 10), Expression (18) is written in the form 








cor. oe n 2 Pp \ sin @, ¢ — sin Xe, ¢ 
Gp Tw ®p! @pt max | (1 ee Kept ne. 
2 

—~Sy (Pmax —1)[1— 


sin @, t — sin X,, t 


d=—‘ja,t : I} . 





Bearing in mind that (15) holds, and substituting the numerical values from (15) into Equation (19, we 
obtain the curve ~¢ (t)/Gp for different values of Prax (Fig. 12). 


It is clear from Fig. 12 that, for Gp < 0, the quantity — € (t)/Gp will be positive. In this case, the trans- 
lent response xy (t) will be aggravated. Its overcompensation, a, , will be approximately equal to the sum of 
the overcompensation of the function x(t) and the maximum value of € (t). Now the graph of oy = f(P,,,,), 
constructed in [1], may be reconstructed taking into account the effect of the lag, r. We thus obtain the 
family Or yy = f(Pmax) for various values of Ge (Fig. 13). For G, = 0, the curve Or yy = f(Pmax) coincides 
with the curve for the ordinary system without lag. We note that for positive values of G at low frequencies, 
the compensation in systems with lag is low compared to the same system without lag. But, in first and second 
order static and astatic systems, the value of G cannot in general be positive. 


Now we turn to the influence of lag on the mag- 
nitude of the regulation time, T. From the transfer 
function x(t) [1, 3] and from ~—e (t)/Gp for the corres- 



































C % ponding typical characteristic (Fig. 10), we obtain the 
P transfer function of the system for various values of the 
ot! a lag, X_ (t). Based on these, the curves T,,, = f(Pnax) 
4 for different values of G, are constructed (Fig. 14). It 
‘ S04 is clear from Figs. 13 and 14 that, for given values of 
30 = 
re ce rg Ory, and T,,,, an increase in |G,| lowers the attain 
7) % able value of P,,,, and increases the required value 
et * of the cutoff frequency we, i.e., the system is subject 
M to more stringent requirements. Thus, Figs. 12-14 
Panini give an approximate idea of the properties of a system 
ou 2 8 with lag elements. 
Fig. 13 On the basis of what has gone before, we can 


assert the following sequence of analytical steps. 
1. We determine w, from the given frequency characteristics, KW(jw), not taking lag into account. 


2, From either the circle diagram or Table 1, and being givenr , we determine the value of G,, which 


has the largest absolute value in the interval 0< w= we (we assume that the point C is the first point of inter- 
section), 


3. Knowing G,, one can immediately estimate from Figs. 12-14 the approximate influence of the lag 
on the fundamental properties of the control system. 


It is known that between the stability margins in modulus 7, and phase y, and the magnitude of the 
compensation there exist definite relationships [1]. By analogy with the usual method, for the cases when 


G # 0 graphs of y, = f(0,) andi, = f(o,)(Figs. 15 and 16) may be constructed for determining the nec- 
essary yr and] , from the condition that o, not exceed a given value. 


We establish from Figs. 15 and 16 that, for a given o Ty’ 2 increase of |G,| requires an increase in 
the stability margins in modulus and in phase. 
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4, Method of Synthesizing Compensating Devices for Systems with Lag Elements 





It is known that, to attain adequate accuracy in reproducing the controlling effect, g(t), which is differ- 
ent from unity, in the frequency interval 0< w< wg, in which the frequency characteristic of the effect is it- 
self sensibly different from zero, the system's frequency characteristic must satisfy the approximate equalities 


{1}: 


P.(o)~1, Q.(@)=0 for OC aay, wg < wy. (20) 


From (7) and (20) we obtain 





1—G.— 
i—G —- 
P (a) = 7 we ” “P for 0 Q@< ay. (21) 
ae ee 
bas 


Expression (21) indicates the desirable form of the function P(w) in the interval 0 =< w= Wg with lag 
present in the system. 


It is clear from (21) that, at low frequencies and 
with sufficiently low values of r or G,,, the function 
P(w) #1. Consequently, the desired fogarithmic gain 
characteristic at low frequencies does not differ from 
the desired logarithmic gain characteristic of the ordi- 
nary system with the identical requirements on the 
accuracy of reproduction of the controlling signal. 


The form of the desired logarithmic gain charac- 
teristic in the intermediate frequency range is deter- 
mined by the basic requirements for stability margins 

*max and from the given desiderata for the system's proper- 

ties for one step-wise controlling section. The desired 

Fig. 14 logarithmic gain characteristic is found from Figs. 12- 
F 16 and Table 1. From the given values of o,,,, T; 
and +, it is possible to determine, from Figs. 12-16, 


the values of w., G_ and the necessary values of stability margins in modulus,! ,, and in phase y,. Itis 
clear that, for the usual linear system, G. = 0. 





fo tw 2 Is 












































db 
x * 
ia” 
~~ 0 Btw 
7] et a Pied I SS 
2” 6,005 - Pe. deen — 
6,7 405 
oe 4, % 
20 JO 40 20 30 « 
Fig. 15 Fig. 16 


In order to satisfy the system requirements for large values of 1, it is necessary that the modulus and phase 
margins be very large, therebycomplicating the construction of a compensating device. 


The form of the logarithmic gain characteristic at high frequencies does not esentially influence the sys- 
tem performance, 
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Fig. 17 


In this manner we determine the desired logarithmic gain characteristic of the system with lag elements. 
After having thus determined the logarithmic gain characteristic, one may easily find the form of the com- 
pensating device, 


Example of the Synthesis of a Compensating Device for a System with Lag 





We assume that we are given an object with transfer function 


1 


KWo (9) = 50s +1) 0.005 +1) O0ls +1) (0.0058 1) (22) 








and that the given conditions on system properties may thus be stated; 
a) the system must have first order astaticism, and a velocity error coefficient C, = 0.005; 


b) the compensation, Or » for one step-wise controlling action must not exceed 40%, and the transient 
response time, T, must not be greater than one second. 


For the given values of control time and compensation, we choose, using Figs. 13 and 14, an approxi- 
mate value of the cutoff frequency, w, = 11.5, Further calculations will be made on the basis of three dif- 
ferent values of r= 0, 0.005 and 0.01 sec. We find the values of gy = rT we for the stated values ofr (Table 2). 































TABLE 2 
t sec. e- deg. Ge. 4,2 db | vz deg. @e @ 
0 0 0 412 35 3 0.19 
0.005 3.3 —0.05 +16 41 2 0.12 
0.01 6.6 —0.1 +19 48 1.2 0.078 


Above all, we must obtain the value of the distortion function |G.| from Table 1, entering the table 
with g and y,. The magnitude of y,, we choose, as is customary, approximately in the interval 20°%<y_= 80° 


In this interval we seek the lowest value of IG] (Table 2). 


Knowing G, and of it is possible, by using Figs. 15 and 16, immediately to determine the required 
stability margins both in phase and in modulus (cf, Table 2). 











Using Formula (10) we find, as for the ordinary linear system, the greatest value of the conjugating fre- 
quency, w:, for which y, =35° (Go. = 0), Yr =41° (G_ =—0.05) and y, = 48 (G_ = ~0.1). 


The values of the mating frequencies ws and wu, are given in Table 2, 
It is easily verified that the stability margin in phase, for w = we, is satisfactory. 


Figure 17 gives the desired logarithmic gain characteristics for r = 0, 0.005 and 0,01 sec. The same 
figure gives the logarithmic gain characteristics for the compensating device at low and intermediate fre- 
quencies, To avoid complications, the high-frequency portion of the compensating device's logarithmic gain 
characteristic was not corrected, which does not essentially influence the system's behavior, 


Thus, we obtain the expressions for the frequency characteristic of the compensating devices: 











for r = 0 
(0.33j@ + 1) (OAj@ + 1) 
KW 4(/®) = “6 26ja 4 1)(0.0055)@ + 4)’ (23) 
for r = 0.005 sec. 
Wo) — —0-ni@ +1) 0-7 + 1) 
KW4 (1) = “8 Baja +1) 0.0055;0 4-1) ’ (24) 
for r = 0.01 sec. 
: (0.83j@ +- 1) (O0.1j@ + 1) 
KW 4 (je) = (12.87@ + 1)(0.0055j@ 4+ 1) ° (25) 


It is clear that the compensating devices for the first and second cases are more easily realized than in 
the third case, 


Appendix 
For greater accuracy in computing ¢ (t) and x, (t), it is necessary to determine the function G. 


Using (4), we put F in the form 





F = B+ Qsinto, (26) 
where 
KW (ja) 
Q= Im T+ KW (je) ” (27) 


U (cos tw — 1) 


B= SUF +V (28) 





with Q defined by means of the imaginary circle diagram, After transformation, (28) may be written as 


1 — cos Tw 


2 i—costo]f,  1— cos tw 
[v+i4 eral 2B = \(2+=3" (29) 


Letting D = (1-cos rw)/2B, we obtain from (29) 











(U0 + (1 + PP + V2 = D(2 + D).z (30) 


Equation (30) is the equation of the family of circles, symmetric with respect to the imaginary axis and 
with their centers on the real axis (Fig. 18). 


If one constructs, on this same plane, the fre- 
quency characteristic of the open-loop system exclu- 
sive of the lag, then the points of intersection of the 
circles with this frequency characteristic give the 
curve D(w). Knowing D(w) and r , one can speedily 
find G( w) using the following relationships: 








(a) 


B 1 — cos Tw 
PF 
F=B+Qsin Ta, C=Tr PF 


D-2 
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SUMMARY 
Fig. 18 


1. The introduction of the distortion function, 

G, allows the carrying out of the property analysis, and 

the synthesis, of compensating devices for automatic 
control systems with lag elements. The idea behind the method of introducing the function G is to establish 
an immediate connection between the frequency characteristic of the open-loop system exclusive of lag, the 
magnitude of the lag, r , and the properties of the system. Knowing the frequency characteristics and the 
value of r , one can, by this method, quickly determine the approximate characteristics of the transient re - 
sponse, X (t), and also estimate, without complicated computations, the influence of the lag, r , on the in- 
dicated control system properties. 


2. The exposition of the method of synthesizing the compensating devices for systems with lag is based 
on the method for synthesizing compensating devices for ordinary linear systems [1]. For systems not possess- 
ing lag elements, we consider G, = 0. The basic circumstance distinguishing systems with lag from those with- 
out is that the curves Trm* f(Pmax), Oru ™ f(Pinax)> Yr = f(a,) andi, = f(o,) are not identical, 
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INCREASING THE SPEED OF CERTAIN AUTOMATIC CONTROL SYSTEMS BY 
MEANS OF NONLINEAR AND COMPUTING DEVICES 


G. M. Ostrovskil 


( Moscow) 


This article considers the questions involved in introducing 
nonlinear devices in certain automatic control systems in order 
to improve the quality of regulation given. 


The basic results in the domain of using nonlinear devices in automatic control systems described by 
n‘th order differential equation, with certain limitations on the variability of one or more coordinates, were 
obtained by A. A. Fel'dbaum [1, 2], A. Ia. Lerner [3] and L. S. Pontriagin, R. V. Gamkrelidze and V. G. Bol- 
tianskii [4]. In the present work, the method investigated in [5] is extended to cover certain automatic con- 
trol systems described by n‘th order differential equations. 


Let the equation for the error x (xX = 6 in ~ @ oy) Of the system have the form: 


2”) + ag) +... +a, 27+ 4,2 = 0 (1) 


We assume that the parameters of the controlling portion vary within limits corresponding to values a;; 


a<a <b, (i= 4,2,...,n—4). (la) 


We choose the parameters aj = bj so that the transient response will not be over-controlled, but will be 
characterized by a long control time. It will be assumed that, in this case, the characteristic equation of 
System (1) will have n real, negative and distinct roots; if the parameters aj = cj, the characteristics of the 
transient response are reversed. 


If, by some nonlinear compensating device or another, we are able to arrange matters so that at the 
beginning of the transient response we have aj = cj, and at the end aj = bj (i = 1, 2,..., m-1), then this will 
assure high speed in attaining the equilibrium state at the beginning of the transient response, and a high de- 
gree of damping at the end. The system's compensating device must work in such a manner that, at a certain 
moment of time, the coefficients aj rapidly change their values. It is clear that if this switch occurs too soon, 
the system will be subject to a long transient response, whereas if it occurs to late, the system will be heavily 
over-controlled. 


One may formulate in the following way the problem which will be solved; at the first stage let the con- 
trol system be described by the differential equation 


2” tea) +... tone t+2=0, 





and at the second stage, by the equation 


ob +... +b, e+ r= 0. (3) 


It is necessary to find the precise moment for transferring from Equation (2) to Equation (3) for which the 
system will not be over-controlled, but for which the control time will be a minimum. It is assumed that dur- 
ing the transition from Equation (2) to Equation (3), the first (n-1) derivatives of the output quantities vary 
smoothly. 


We shall consider those systems which, for a unit excitation, have the following initial conditions; 


ty <0, 2) >0,...2%" >0 (4a) 


or 
i, as a (n—1) 
Lo P* 0, Lo _— 0, eee Lo — 0. (4b) 


If the system has initial conditions of the form (4a), then a necessary and sufficient condition for mono- 
tonicity of the integral curve will be 


F (Gq - - - Zn) = Yaa - + - YnTi + (Yas - ~~ Yn—1 + Yes +++ Yn—sY¥n +--- 
vee Hb Ys¥a- ++ Yn) Le + (Yas - - - Yn—2 + Ya¥s + - - Yn—s¥n—1 + Yo2¥s - - - 
+0 + Yn—a¥n—2¥n—1 +--+ + Ya¥s- ++ Yn) Zp... + (yaya + Yo¥a + - 
oes EF Yan + Ya¥a t+ +++ + Ya¥n bt - - s + Yn—2¥n—1 + Yn—2¥n + 
+ Yn—1¥n) En—2 + (Yo + Ys +--+ + Yn) Zn + In <0, (5a) 


where x; = X, X, = Xt scoces Xp = Soa (Cf. Appendix). 


In the sequel it will be necessary for us to consider the hyperplane 


F (2, %q, . +» %n) = 0. (5b) 


It will be assumed that, in x;, X2,.....X, phase space, the initial points, determined by Conditions (4a), 
lie “beneath” Hyperplane (5b). We shall say that a point with coordinates X;, X2,....Xp lies “above” a surface, 
Kn = £(Xq,Xg,-00yXp-g), If W (Ky, Xp,....Xn-1) = Xp — f(Xy,Xe,....Xp-) > 0, and lies “below” this surface if 
¥ (Ky. Xg,-.-5% 1) < 0. 


The desired system will be such that a; = cj so long as the phase point lies “below” Hyperplane (5b), 
and aj = b, for all points lying on or above this hyperplane. Mathematically, this is written: 


c;, F (a, z,...,2° "<0, 


hea 7 : (n—1) 6 
Oe F (6. 2...2,8 .. >>. (6) 


We now consider the transient response with initial Conditions (4a). 
n-1) 


Since, as was assumed, F (x9, Kgreee Xt 
following differential equations; 


< 0, then aj = cj and the control system is described by the 
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We shall assume that the portion of the hyperplane 


C1En + Cop +... + 4, =O, (8) 


lying in the region D(x, < 0, X, =0,...,x,, = 0) lies “above” the portion of Hyperplane (5b) lying in the same 
region, This can be attained by making cj sufficiently small. We shall prove that, during the switching pro- 
cess, the following inequalities will hold; 


z,<0,2,>0,..., & >0O. (9) 


In all cases these inequalities will hold until the phase point falls “below” Hyperplane (8) since, in this 
case, X, > 0. 


To the extent that Hyperplane (5b) lies “below” Hyperplane (8), the sooner will the integral curve in- 
tersect Hyperplane (5b), this intersection occuring in region D, It is easily seen that the integral curve cannot 
enter the region x; > 0, being “below” Hyperplane (5b). 


Indeed, while the phase point is “below” Hyperplane (5b), the inequality, F (x,,Xg,....X,-3) < 0, must 
hold. For, if the integral curve intersected the hyperplane x, = 0 then, at this point, we would have x, = 0, 
Xg > 0,...,X, > 0, from whence F (xy,Xg,...,Xn) > 0, contradicting the above. Thus, the system cannot be 
over-controlled at the intersection of the integral curve with Hyperplane (5b). 


After switching, the coefficients a; will become equal to bj, and the system will be heavily damped. 
Hence, the system, while tending towards an equilibrium state, will not be over-controlled (see Appendix), 


If switching occurs too late, the system will certainly be over-controlled and, if too early, the system 
will be heavily damped during a long interval of time, thus increasing the duration of the transient response, 
In the case Xg> 0, X9< 0,....x(272) < 0, a curve whose initial point lay “above” Hyperplane (5b) would never 
leave the semiplane x; > 0, i.e., F (xy,Xg,....Xp) = 0. In order to obtain a device which would function cor- 
rectly for both xg< 0 and Xg> 0, it is necessary to take the sign of x into consideration. Thus, a control system 
will now be described by the following equation: 


og” +a" +... 4+a,.7+2=0, 
where 


sa bey, SE tte Fen rn Or” KO (10) 


. . i=i, as gad —_ . 
7 b;, 2F (z, z,... 2") >0 ¢ =< 


This means that, in an automatic control sh ei there must be included an analogue computing device 
which would generate the function xF (x,X,...,% ). 

















The schematic for such a system is given in Fig. 
1, At the input of the computing device, 1, there is 
impressed the error signal, x, and at its output appears 

a 6 the f ~ (n=1) * 

vn "Mee out e function x F (x,x,...,x ), controlling clutch 3; 
by means of a linear compensating device, 2, it is pos- 
sible to change certain of the coefficients, a,; the con- 
trolled object, 4, is also included in the system. 


























If clutch 3 is open, the subsidiary feedback path 
will not exist, and the system is described by Equation 
(2). If clutch 3 is closed, then certain aj will become 
Fig, 1 equal to bj. Opening and closing of clutch 3 will de- 
pend on the sign of the function xF(x,X,....x(?)), 











It is clear that the output of the computing device may be fed to several clutches, controlling the work 
of several compensating devices. 


We have cited a class of systems which, for unit initial excitation, have initial conditions of types (4a) 
and (4b). If the transfer function of the open system has the form: W(p) = k/Q(p), where Q (p) is a poly- 
nomial, then the initial conditions will be xg =—1, Xgq = X39 = ....Xpo9= 0. If the open-loop transfer function 
has the form W(p) = (Tp + 1)/Q(p), then the initial conditions will be xg=—1, X99 = ...= X(n-1)9 = 9 Xn> 0 [6]. 


The latter type of system underlies many systems of isodrome control (i.e., strict control of rotational 
speed under varying load conditions), used in the regulation of production processes, 


It is easily seen that, in the first case, during an infinitesimally short interval of time, x, > 0,x,_, > 0, 
soeXg > 0, X1 < 0. Indeed, integrating System (7) with the initial conditions xy9 = —1, Xg9 =....=Xp9 = 0, over 
a sufficiently small period of time, we obtain 


%, =t>0, 
12 

2-1 = a > 9, 
3 

Zan. = 3I > 0, 


where t is a small interval of time. 


Reasoning further as in the proof of Inequality (9), we may show that the phase point must intersect 
Hyperplane (5b) and, once upon it, tends to an equilibrium state, The same reasoning is applicable to the 
second case. It is only necessary that the phase point, 
determiined by the initial conditions lie “below” Hy- 
perplane (5b), which can be attained by properly choosing 

PU serpy out the quantities y,....’n. It should be mentioned that 
i & e—— s the applications of the method given are not restricted 
to these two classes of automatic control systems, or 
to systems characterized by initial Conditions (4a) and 
u | (4b). In each specific case when Conditions (4a) and 


(4b) do not hold, it is necessary to investigate the con- 
trol system. 
































Fig. 2 We consider, by way of an example, a third order 
control system (Fig. 2), The invariant portion of sys- 
tem 1 consists of the regulated object and the executive 

mechanism. The transfer function of the open- looped system has the form: 
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Ww an 
(”) = Fp +) Tt): (11) 





For unit excitation of the given system, the initial conditions will be x, = —1, X—= 0, X9= 0. 


We denote by n the feedback loop's speed coefficient, and by m* its acceleration coefficient. In the 
case where there are no nonlinear compensating devices, consisting of clutches 2 and 3, in the link consisting 
of the subsidiary feedback loop and computing device 4, the system equation will be 


[7,747 yp? + (7,7 4+ 731 +m p+ (7 +n) p+ KOg., = HO, (12a) 


If we introduce the nondimensional time 


Vix 
t=) fT." 





then Equation (12a) takes the form: 


Gur + 48 our + Bout + ur = Pin (12b) 


where 





Ps TiTs3 ot T.T3 + m? B= Ts +n 


3 


A 
3 
VY kT? 27? VT ,T.T; (13) 


and the dots represent differentiation with respect to r. 


For m = n = 0, the coefficients A and B take the values; 


i T:T3 a T2T3 Ts 


Ao ’ B= 3 


; ; 
meena a 13a 
V ererer? VET TT, o— 


For some fixed m and n, the coefficients A and B will take the values Am and Bp. 


The nonlinear compensating device will function in the following manner. Computing device 4 gen- 
erates the function u = © (x,x,X'), where x = @jn ~ @oyt- The signal u is fed to clutches 2 and 3, which pass 
the signal through the subsidiary feedback loop if u =0, and inhibit itifu < 0. 


The control system will be described by the following equations; 


@+ Agr + Bor+2=O0 for u<0, (14a) 
2+A,2+ B,2+2=0 for u>0. (14b) 


Let the roots of the characteristic equation be — a, —y +iw in case (14a), and ~y; ~y,-7 s in case (14b), 
with y1<y2<7s3- 











The function at the output of the computing device will have the form; 
u = @ (x, 2, 2) = 2 [ysyer + (x2 + ys) 2 + 2]. 


Equation (14a) has the following solution: 





w(t) = Cye** + Cyre—* sin wt + Cye~ ™* cos wt, (15) 
where 
y+ o? _ (a? — y? + ay) _ _4(a— 2y) 
A Gaya = oa Fal] = Gat ah 


and Equation (14b) has the solution: 
z(t) = Bye~™* 4+ Bye~™* 4+ Bye**, (16) 


where 








B, = 12 + (y¥2 + Ys) %1 + 2 By = Ut + (1 + Ys) 21 + 2; 
(y2 — Y1) (ys — Y:) ; (v1 — Ys) (ys — Y2) 
By = 1%" +r + ye) 21+ 2 ; (17) 
(Y¥1 — Ys) (y2 — Ys) 





where x, ,%,,%; are the values of the coordinates and their first and second derivatives at the moment of switch- 
ing. 


If switching occurs on the surface @(x,x,X) = 0, in this case B, = 0, since the numerator of B, coincides 
with the quantity in brackets in the expression for ®(x,x,%) for x = x,, X= X, and X = X. 


We find the expression for u by substituting the value of x from (15) in @ (x,x,X). 


u = (Lye~** + Loe sin wt + Lye coswt)z, 


where (18) 


Ly = C; [ysvs — 4 (v2 + Ys) + @?], 
Ly = ya¥sCa— (Cay + Cee) (2 + ys) + 2ywCs + C2 (y* — w?), 
Lg = y2¥sCs +- (2 + Ys) aC; — aC. 


Let us assign the numerical values; T, = 35 sec., Tz = 6.55 sec., T, = 3 sec, and k = 0.5. Then Ag= 2,02 
and By = 0,54, We set x= —1. 


The solution (15) has then the following form 


a (t) = —0,112e—** — 0.319e—* sin 0.7077 — 0.888¢—"* cos 0.707. 


Let the speed and acceleration coefficients of the feedback loop equal; n= 4.55 and m = 15.8. Then 
Am = 6, Bp = 8.5, Hence, y; = 0.13, y2 = 2.0, y3 = 3.87, Ls = — 7.74, Ly = 1.4, Ly = 0; the function u vanishes 
at time 7, = 1.97. Hence, x (r,) = — 0.46, X(r,) = 0.57 and X¥(r4) = 0.21. From (17) we get that B, = 0, 
B, = 0.65 and B, = 0.19. 
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The time, rg, during the course of which X(1) becomes less, in absolute value, than 0,06, will equal 
1.17. Consequently, the time for the transient response equals T = rT, +7, = 3.14. 


In this case, the curve for the transient response 
(Fig. 3) is DEyG,. On are DE, of the curve we have 
the solution to Equation (14a) and, on arc E,G,, the 
Z(t) f Z J 4 g U t solution of Equation (14b). 





For comparison, we take the curve DEG, of the 
transient response of a third order linear system in which 
A = B = 3 and, consequently, y; = yg = ¥3 = 1; the equa- 

: 











tion of the curve DE Gg is: x(r) = 1 +T + - :™ 


Fie. 3 In this case, the time for] x( 7 )| to become less than 
8: 0.06 will equal 6. Thus, in the nonlinear system, the 
time of the transient response is almost twice as low. 


We sought the moment of switching under the assumption that qj and bj were given, To give a rule 
which would unambiguously determine the coefficients c; and bj is difficult. However, it is possible to show 
that, in the first period, the smaller the coefficient cj, the quicker does the system move. The coefficient by 
may be found by constructing, in the aperiodic region, curves of equal stability, using the method of Tsypkin 
and Bromberg, and choosing bj by the curve corresponding to the least degree of stability. 


Appendix 
We consider the equation 


2) +. dyel™—D 4 bol) 4 4b, ete = 0. (A.1) 


Assuming that the characteristic equation has n real, negative, distinct roots, —y ,>—y >....> “Yn we 


find the necessary and sufficient condition that the integral curve of Equation (A.1) be monotonic for the initial 
conditions xX» < 0, X—> 0, X»> Ovex 2” 9 > 0. We turn to the system of equations; 





dz, 
yo — bx, — bye, — ++ — By 
dz, _, 
= £ . 
at bd (A -2) 
dz, 
a = 2. 


We consider x; ,Xg,...,X,, phase space. In it there are, among the integral curves of System (A.2), the n 
lines; 


—yt 
m14=>——e 
7,8 


—yt 
2,=—yje v4 ’ 


(A.3) 


= (—1)"yi"—"e “' (i = 4, 2,..., a). 


We conceive of the lines (A.3) as the axes of a new oblique coordinate system, Uy, Ug,....Up,, where we 
consider, as the u, axis, the line corresponding to y;, as the ug axis the line corresponding to yg, etc. As the 
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positive semi-axis of the u, axis, we take the semi-axis lying in the half-plane x, < 0. The formulae for 
passing from the uy ,Up,...,Up system to the x;,X3,...,X, system have the form [7]; 


z= —(u,+u,+.-..+4,), 
By = ¥y% + gM", +--+ YnY ns 
Zs = —(¥j4, + YeMgt--- + Ru), 
3 1 2°2 nn (A.4) 


——-: @ = Gri « « Bee TF eye « ws 


a, = (—1)" (yf *u, + ypu, +--+ yaw). 


In what follows, we shall need the expressions for u, and ug in terms of x;,Xg,...X,. Solving System (A.4), 
we get 





es (Y%; - ¥;) 
penn eran (¥; —¥;) [¥g¥g- ++ Yuta + (Va¥g- ++ Yn—1 t Ye¥a- ++ Yn—2¥n + 
1<j<ci<n : . 


+ Yes + ** Yn—s¥n—1¥n + + + > + eV: -* Yn) 2g + +++ + (Ve¥3 + - Yn—t4 + 
+ Vols ++ + Yn-i¥n—~i¢e T+ ++ Ye¥ar + * Yn-ite F 
+ ¥a¥q ++ + Yn—t41 Yn—ita H+ + Vals +++ Vn—ing H+ F Vig Vien ++ Yn) i + 
+ +t (gig t+ Va¥g + +--+ Ya¥n + Ys¥q +--+ ++ Ys¥n +--+ + Yn—1yn) Zn—2 + 
+ (¥g + gy + --- + Yn) Zn-1 + Fnl- 


(A.5) 


The coefficient of xj takes the form of a sum of products of (n-1) of the n roots,|y; | For ug we obtain 


II (%— Y;) 


1<j<icn(t +2; j#2) 


7 (y4— Y;) 
1<j<i<cn 
+ V178"**Yn—3ln—iln t*** + Y¥8¥ar* Yq) Fa Hest (Yt¥8* + niga + V1%9"**Yn—i¥n—it+2F (4.6) 
Hees Hb ¥a¥e*** Yn—ite + Ys¥e °° * Yn—ipaYn—its HF °°* H M¥s *** Yn-its + 
Het FV rVi4g *** Vn) 2% e+ + UV2¥s + Ya¥e + ++ Y2¥n + t8¥e t +++ + ¥8%n + 
tert Yn—1Yn) Tn—~2 + (Ye+ Yst**: t+Yn)*n—-1 + z,,). 





[¥1¥s***Y¥n% + (¥1¥8" + -Yn—y + Y1¥2"** Yanan t+ 


The uy,Ug,...,U, coordinate system will be normal [8], and in it the equations of the integral curves 
will have the form: 


(A.7) 


If uy = 0 at the outset of the motion then, for all time, uy = 0, i.e., throughout the motion the integral 
curve will lie in the hyperplane u, = 0. We designate this hyperplane by the letter,L. We might make a 
similar statement for ug, etc. It is easily seen that an integra) curve which does not lie in any of the hyper- 
planes u, = 0, ug = 0,...,u,, = 0, may intersect them only at the origin of coordinates. 
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Consider hyperplane L. We will show that, if the integral curve lies “above * hyperplane L, it will leave 
the half-space, x, < 0. If the integral curve lies “below” L and, at the initial moment, x; > 0, Xg > 0,...,X, >0, 


then it will not leave the half-space x,< 0. We get the equation of hyperplane L from (A.5) by setting u, equal 
to zero. 


It is known that all integral curves of System (A.2), excluding those for which uw = 0, are tangent to the 
uy axis [8], being tangent to the positive semi~axis for u,g > 0, and to the negative semi-axis for ug < 0. 


Consider the integral curves whose initial points lie in region D "“above* hyperplane L. Substituting the 
coordinate values for these points in the right side of Equation (A.5), we get uyg< 0. From (A.7) we then obtain 
that u, < 0 during the entire motion. Consequently, the integral curve will be tangent to the negative half of 
the u, axis, which lies in the half-space, x; > 0. The integral curve must of necessity then leave half-space 
X, < 0. We have thus proved the first part of our assertion. We now assume that the initial point lies in the 
region D “below"™ hyperplane L. Substituting the coordinates of this point in the expression for u,, we find that 
Uyg > 0, and that, therefore, u, > 0 throughout the motion. Consider the intersection of hyperplanes u, = 0 and 
ug = 0. For this purpose, we rewrite the equations for these surfaces in the following manner; 


F (xy, 22, +--+ 5%) = Ye¥a** “Yq tit [Ye¥s" + Y_ig te +H ysYer + Yn) 2 + 
+ +++ + lysva + v2¥e+ +++ + Y2¥n + (Ysa + +++ + Yp_i¥p)] Zag + (A.8) 
+ [ye + (ys + +++ Y_)i p14 + 2%, = 9, 


Y1¥9 °** Yn%1 + [yr¥s ++ > Ynig + ++ + Y¥8¥a*** Yn) 22 + +++ + Eyre + yrva + + 
2 Eb VIYg + (Ys¥e + ° ++ + in—1¥n)) tng + [1 + (Ya + +++ + Yq) Zp + %_ = O- (A .9) 


In the coefficients of xj in Equations (A.8) and (A.9), identical expressions are isolated within parentheses. 
Outside the parentheses are terms which differ only in that they have y, in the second equation where they have 
Yz2 in the first. We multiply (A.8) by y, and (A.9) by y 2, and subtract the first from the second. With this, the 
terms outside the parentheses drop out, and the terms inside the parentheses may be written as; 


Ys" **Yn(¥2— 1) Za + +++ + ysve + +++ Yn—¥n) (¥2 — Yt) Spe + 
+s ++ bY pM Y2—Y1)Fp_s +(e — Y1) Zn = O- 


We divide through by the common factor (y2,~ 7) to obtain 


Ya °° Yn%at ++ + rsva + +++ + Yn—1¥n) Zn—2 + (18 + +++ + Yn) Fay + Zp, = 0. 


Since all the coefficients of this equation are positive, there is no solution for xg > 0,...,X, > 0. Hence, 
the hyperplanes u, = 0 and ug = 0 do not intersect in region D. 


We now show that hyperplane uy, = 0 is “above” hyperplane ug = 0 in region D. Consider the line parallel 
to the x, axis; 


a= —4, y= 23 = +++ =F = 0, (A,10) 


The points of intersection of this line with Hyperplanes (A.8) and (A.9) will be the points B, (-a,0,...,0,X%n) 
and By (-a,0,...,0,Xn), where Xp = Y 273-7 n4 and Xp = 7173--sJna. We have X_> Xp, since yg > yy. Hence, 
point B; lies above point B, on Line (A.10). Consequently, insofar as hyperplanes u, = 0 and ug = 0 do not in- 
tersect, hyperplane u, = 0 will be “above” hyperplane u, = 0 everywhere in region D. Consider now a point 
of hyperplane u, = 0 inregion D. It lies “above” hyperplane ug = 0, From (A.6) we then obtain for this point 
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that ug > 0. We now show that, if the initial point lies in the region D “below” hyperplane L, then the inte- 
gral curve starting at it will not leave the half-space x, < 0. We prove this by the use of mathematical in- 
duction. It is proved for equations of the second order in [5], hence we may assume that the assertion is also 
true for equations of order n 


In this case, the initial point lies in the region x; < 0, x, > 0,...,x,) > 0 “below” the hyperplane uy, = 0. 
This region in uy,Ug,...,Up, is given by the following inequalities; 


4 + Ug+---+u,>0, 
Yili + Yota + +++ + ¥,4,> 0, 
Vy + Ygue + +++ yyy <0, 


(— 1)? (yf Pan + yg ue + +++ + yA Pu) > 0, 
(— 1)" (yf + yg tue + ++ + ye) > 0, 
u,> 0. 


(A.11) 


An integral curve whose initial point lies in this region does not, according to the assumption, leave 


the half-space x, < 0, 1.e., does not intersect the plane uy, + ug+... +u,, = 0. 


We now consider.a system of the (n+ 1)'st order. The region D is defined by the inequalities 


uy + uz + PAS +u,4,>9, 
yita + Yota + sites + Yn4i4n41> 9: 


ae eS eS o's 6 82 + 8a es 


(A) yes vg a + + RF Mngs) > 9, 
(—1)" 4 (yf + vga + -++ + yh Unis) > 0. 


(A.12) 


Consider a point on the hyperplane u, = 0 in region D. The totality of such points will constitute an 
n-dimensional subspace. These points, lying in region D on the hyperplane u, = 0, will satisfy the inequalities 


Ug + Ug + +++ +41, > 0, 
Yalta + yatta tees + Yn i4nia > 9, 


(— 1)" (yg 7a + +++ + YT Mn 4s) > 0. (A.13) 


We previously showed that, for points lying on hyperplane L, ug > 0. Thus, in the n-dimensional sub- 
space, u, = 0, the points lying in D satisfy the same inequalities as in (A.11), except that now, instead of the 


u, axis we have the ug axis, instead of the ug axis we have the u, axis, etc. On the basis of our previous assump- 
tion, an integral curve starting at one of these points could not intersect the hyperplane ug+ Us +...U, 4 , = 0, 
i.e., could not intersect the hyperplane x, = 0, since, for u, = 0, X; = ~(ug+ Ug +...+U,) 4 4). We have thus 

shown that an integral curve lying on the hyperplane u, = 0 and in region D can not leave the half-space x; < 0. 


It can be shown that integral curves lying in region D “below” hyperplane L will not leave the half-space x, < 0. 


Consider the point A(x;,Xg,...X,) lying on hyperplane L and point A (X;,%g.--.,X%) lying on the line parallel to 
the x, axis and passing through A, where A lies below A. From this, x = Xj ,Xq = Xg,..4.Xp-y = Xp-ys Xp > Eye 
Consider the integral curves starting at, respectively, points A and A. The initial velocity, acceleration, third 
derivatives, etc. of these integral curves are directed towards the initial excitation. Hence, all initial deriva- 
tives up to the (n-2)'nd order are equal for the two curves, but the initial (n-1)'st derivative of the curve start- 
ing at point A is less than the corresponding quantity for the curve starting from point A. However, we have al- 
ready showed that the integral curve starting from point A cannot leave the half-space x, < 0. Consequently, 
the integral curve starting from point A can certainly not leave this half-space. 
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CONCERNING A POSSIBILITY OF DETERMINING AN AXIAL HYDRODYNAMIC 
FORCE IN A VALVE 


I. M. Krassov, and B. G. Turbin 


(Moscow) 


The paper deals with an axial hydrodynamic force that appears 
in valve hydraulic amplifier when the working liquid flows through 
it. The amplifier described may be a meter of the force. The ex- 
perimental results are presented. 


One of the specific characteristics of valve hydraulic amplifiers, which find wide application in auto- 
matic control systems, such as the following and others is the appearance of axial hydrodynamic forces which 
act on the valve from the direction of flow of the work fluid [1, 2]. 


The axial hydrodynamic force depends on the flow parameters of the work fluid (the drop in pressure py 
of the fluid on the working edges of the valve, the fluid density, the valve opening x) and on the constructional 
parameters of the amplifier (width of the through section of the valve, radial clearance between the shoulders 
of the valve and its sleeve, radius of curvature of the working edges of the shoulders, axial distance in the valve 
between the centers of the inflowing and outflowing stream, etc.). This force consists of a steady and an un- 
settled component and, arising only during the flow of fluid through the valve's opening, endeavours to close 
it, i.e., counteracts the force of the control device. 


The effect of the axial hydrodynamic force on the work of the hydraulic amplifier is particularly per- 
ceptible in amplifiers having a large power amplification factor [3], such as those found in electrohydraulic 
servomechanisms. The output power of these is raised by increasing the pressure and consumption of the work 
liquid while the sensitivity is increased by utilizing low power electrical control devices. In this case the 
appearance of the axial hydrodynamic force imposes much more severe requirements on the electrical control 
devices. Thus, a study of the nature, magnitude, and effect of the axial hydrodynamic force on the operation 
of the hydraulic amplifier, as well as a study of methods to minimize it is of definite interest. 


In the course of testing one type of two-stage hydraulic amplifier with a high power amplification factor, 
the feasibility of determining the magnitude of the axial hydrodynamic force from the readings of a mano- 
meter on the needle chamber was exhibited.. This possibility results from the principal characteristics of the 
amplifier itself, a diagram of which appears in Fig. 1. In this amplifier the displacement of the needle 4 brings 
about a change in pressure in the pin chamber A. The difference between the pressure of the fluid on the face 
of the valve 1 in chamber A and the pressure of the fluid on the face of this valve in chamber B together with 
the force of spring 5 causes the valve to be displaced. The fluid is able to enter one or another working cham- 
ber of the output mechanism, bringing the piston of the latter into motion. When the flow commences, an 
axial hydrodynamic force arises on the valve and persists until the end of flow of fluid through the opening of 
the valve, i.e., until the end of motion of the output mechanism which is controlled by the amplifier. 


We will assume that the forces acting on valve 1 are positive when they act to the left (from chamber 
A to chamber B); we will also assume that the opening of the valve x is positive when it is displaced to the 
left of its neutral position (towards chamber B). 
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Fig. 1. Principal diagram of the hydraulic amplifier: 1) valve; 2) 
sleeve; 3) throttle; 4) needle; 5) spring. A) needle chamber, B) 
sink chamber. 


Under these conditions the equilibrium equation of the valve, in the case of an established flow of fluid through 
it, has the form: 


(N +cz) + Fp +R=F,, (1) 


where N is the initial tension of spring 5, corresponding to the neutral position of the valve, c is the stiffness 
of the spring, x is the opening, Fa, Fp are the work fluid forceson the faces of the valve in chambers A and B, 
R is the steady component of the axial hydrodynamic force. 


Force Fg may be assumed to be constant, whereas force F, is defined by the equation 


Fa=kp, (2) 


us 





where p is the pressure in the chamber A, k = (D? — d?) is the constructional constant of the valve (D Is 
the diameter of the valve's shoulders, d is the diameter of the needle opening in the valve). 


After the cessation of flow (caused, for example, by the stopping of the output mechanism) force R will 
disappear, equilibrium of the valve will be disturbed and, the position of the needle remaining unchanged, 
the valve will be displaced further. The pressure in chamber A will change; the equilibrium equation will take 
on the form; 


(N + cz’) + F, = F,, (3) 


where x' is the new opening of the valve and Fa", the force due to the liquid pressure on its right face. 


Subtracting Equation (3) from Equation (1) and utilizing (2), one can obtain the axial hydrodynamic 
force R: 


R=k(p—p’)+c(z#’— 2). (4) 
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The difference (x‘'~ x) can be determined through the pressures p and p’ with the help of the static charac- 
teristic of the first stage of amplification, which is determined experimentally. This characteristic represents 
the dependence of the pressure in the needle chamber on the displacement of the needle with respect to the valve 
and, in this particular case, is well approximated by the equation 


vat Pst 
P= 741 ° (5) 


where por is the static pressure in needle chamber A; p, the displacement of the needle with respect to the 
valve, measured from the position at which the needle completely fills the opening into the valve. 


Taking into account that 


ip’, kg/cm" 


d x’ —2 =p’ —p, (6) 


where p * is the displacement of the needle in the ab- 
sence of flow of work fluid through the valve, p is the 
the same but during the flow of work fluid through the 
valve, we will obtain utilizing (5) the following ex- 

pression for determining the magnitude of x'— x: 





I fog Nya 


Fig. 2. Pressure in needle chamber as a 
function of control current. p) during 


and p*) in absence of flow of work fluid 2’ —z=0,134 (\/ 2e-s Pst __ —\/ *—1) — 4 (7) 
through valve. ' 


It follows that by measuring p and p‘ in the needle chamber and knowing the suffness of the spring c 
and the constructional constant k, it is possible to determine from (4) and (7) the magnitude of the axial hy- 
drodynamic force acting on the valve together with the frictional forces. 


Figure 2 shows the results obtained from measurements of the pressure in the needle chamber in the pre- 
sence and absence of a flow of fluid through the valve, at different control signals i causing a negative dis- 
placement of the valve. Kerosene was used as the work fluid in the experiments. The supply pressure on the 
entrance to the second stage of the amplifier (the valve) was equal to 50 kg/cm* The experiment was carried 
out for a valve having shoulders with a diameter D = 25 mm, needle opening diameter d = 2.5 mm, radial 
clearance equal to 0,006 mm, The stiffness of the spring 5 (Fig. 1) was equal to 3,45 kg/mm. 


The experimental and calculated data necessary to determine the axial hydrodynamic force are tabulated 
below: 























i,ma | p’,kg/cm*| p. kg/cm’| *(p’—p). kg | o(x’—2). kg R, kg 
3 2.7 2.5 0.986 0.059 1.045 
5 2.85 2.2 3.204 0.172 3.376 
10 2.8 1.5 6.409 0.455 6.864 
15 2.55 1.25 6.409 0,552 6.961 
20 2.5 1.15 | 6.655 0.621 7.276 


Figure 3 shows the relationship between the axial hydrodynamic force and the valve opening. The distri- 
bution of the experimental points indicates that the accuracy of the experiment was insufficient. On the strength 
of this the given relationship must be considered as approximate, From our point of view, the important fact 
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Fig. 3. Axial hydrodynamic force as a 
function of valve opening. 


is that a two stage hydraulic amplifier of this type 
(needle-valve) may serve as a meter of the axial hy- 
drodynamic force. Such a method for determining 
the axial hydrodynamic force eliminates all forms of 
loading arrangements. 
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THE THEORY OF THE STRUCTURE OF COMBINER MECHANISMS 


S. M. Iakovlev 


(Moscow) 


INTRODUCTION 


Mechanisms which transform spatial combinations of displacements of input elements (levers, keys etc.) 
into displacements of output elements are widely used in telegraphy, telemechanics, centralized railroad con- 
trol, the textile industry etc. Mechanical coders and decoders are examples in telegraphy and telemechanics, 
signal control panels in centralized railroad control (CRC), and automatic pattern-producing devices for weav- 
ing machines in the textile industry. 


We assign these devices the general name of combiners. Although they are widely used, no systematic 
theory of their structure has so far been presented, which has made analysis and synthesis very difficult, Empir- 
ical methods have been used almost exclusively in design. 


In telegraphic equipment, for instance, a set of gaps and projections is chosen for each key, and the de- 
sign then checked for correct operation and absence of false signals. In signals equipment for centralized use 
the levers and interlocks are mainly developed in the same fashion. The numerous combinations of possible 
routes make the choice exceptionally difficult, The situation is no better in the textile industry. 


The first attempts to rationalize analysis and synthesis have been made in CRC. A routing theory has 
been developed [1] which is termed mathematical because it represents the operation of combiners by symbolic 
expressions analogous to fractions. However, the structure of combiners is not considered systematically, so 
synthesis problems for CRC are still mainly solved by selection. 


In 1948-9 the structure of combiners was studied systematically, and analysis and synthesis were based on 
Boolean algebra (algebra of logic) [2].* Here we give the same results, somewhat extended.* * This theory 
raises some new problems demanding further study. : 


Structure of Two-Position Combiners 





Any combiner, of no matter what design or function, has input (command) and output (effector) moving 
mechanical elements; input displacements are transformed to output ones, In telegraphic transmitters the in- 
put elements are key levers and the output keying mechanisms (cams), while in the decoders the input elements 
are selection racks or rings, and the output letter levers. In train controls the inputs are the signal levers which 
close routes, and the outputs the routes left open. 


In the simplest case each input and output element has only two positions of rest. The mechanism is 
then a two-position one. The two positions are referred to as the normal and the set. 


Some such mechanisms, e.g., the ST-35 transmitter in telegraphy, transform the displacements of the 
input elements directly into displacements of the outputs. The keying racks in the ST-35 transmitter are fitted 
with equilateral or right-angled triangular teeth on the key panel side, for instance (Fig. 1). When the push- 
button P is pressed the bottom of the lever presses on the racks and causes them to move in directions deter- 
mined by the slopes of the sides of the triangles. 

*The results were presented at a symposium on mathematical logic at Moscow State University in 1949. 
* *Certain others have also applied Boolean algebra to combiners (e.g., [3 and 7]), but in another way. 
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In other such mechanisms (e.g., in the Creed telegraphy transmitter) the input keys only prepare the out- 
put elements for displacement (free or lock them), while the actual displacements are effected by other com- 
ponents, 


For instance, the press-keys P in the Creed (Fig. 2) lock the rods R and stop their displacement by the 
springs S when the reciprocating head H frees the rods, In railroad control panels some route levers block 
some displacements of the others. 


Combiners of the first type differ from those of 
second type only in design. Any one of type one can 
be made to operate in two steps (preparatory, when 
the inputs are displaced to free or lock the outputs, 
and working, when the outputs are displaced by some 
drive mechanism) by changes in detailed construction 
only. This then gives us a type two mechanism of the 
same structure as the original mechanism. 


In future we shall therefore consider only type 
two, in which input displacements merely prepare 
output ones, 





seatial In the traditional notation [1] a combiner is 


denoted by a kinematic system in which the input, 
intermediate and output elements are represented by horizonta! and vertical lines (levers) hinged together 
via triangle links. Where coupled levers intersect a lock is shown (i.e., a part which blocks the movement). 
Thus Fig. 3 shows that element A locks element B, and Fig. 4 that element B locks element A. This convenient 
symbolism will be used in future. 





Fig. 2 


Algebraic Form of the Locking Conditions 





We use Boolean algebra to study two-position combiners; this algebra has been used successfully with 
two-position (relay) electric circuits [4] in structure studies,* 


Let A, B, C,... be elements of the combiner. Let a, b, c... be variables representing their positions 
which take two values (0 and 1) only. Let the zero correspond only to the normal state, and unity only to 
the set state. 





*The fundamentals of Boolean algebra are presented in ch. 1 of [4], and also (using other symbols) in (5) in 
the chapter on the calculus of prediction; also in [3]. 
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Then the condition that some element X be locked in its normal position by A, B.C... can be expressed 
as the Boolean function f(a,b,c...) in the variables a,b,c.... 


Combinations of a,b,c,,.. correspond to combinations of A,B,C.... Suppose f(a,b,c...) = 0 if the com- 
binations correspond to combinations which lock X in its normal position, and f(a,b,c ...) = 1 if the combina- 
tions correspond to combinations which free X for displacement to the set position, 


Then the condition that some element X be locked in the set position can be expressed as the Boolean 
function g(a,b,c...) which is zero if the combinations correspond to combinations which lock X in the set posi- 


tion, and unity if the combinations correspond to combinations which free X for displacement. 


In this type of mechanism the elements will be free in both positions for the same combinations of the 
locking elements, In other words, displacement of one element does not displace the others and does not alter 
the interlocks which control itself, so it can return to its original position. Hence f(a,b,c ...) = 1 when and 
only when g(a,b,c ...)= 1. 


Not all possible combinations of positions actu- 
ally occur when the mechanism operates. Unused states 
occur as in relay systems [5], i.e., unused combinations 

ie of positions. The locking and locked elements exchange 

4 = roles in real combiners. For example, while lever A in 
the set position locks B in the normal position, lever B 

¥ od in the set position locks A in the normal position. 


f and g are undefined for the unused combina - 

4 686 A 86 tions i.e., can be chosen at will. Choice of values of 
the unused combinations only affects the external form 
of the functions and has no effect on the values for the 
used combinations, which are the only ones which de- 
scribe the real locking conditions, The condition f(a,b,c,...)/, - ky, b = kg, ¢ = ky,... = % where ky, kg, ks,... 
are constants(i.e., zero and one) and “/% is the substitution symbol, implies that g(a,b,c,...)/a = k,, b = kg, 
c = ksg,... is undefined, Conversely, the condition g(a,b,c,...)/, = ky, b = kg, c = ky,... = 0 implies that f(a,b, 
Cyeee)/y = ky, b= ke, c= Ne seco is undefined. 




















Fig. 3 Fig. 4 


Hence f = g except as regards the unused states. If we put that f = 0 when and only when g = 0, we have 
f = g precisely. 


This condition on f and g for the unused states is not necessary, but is convenient. 


Algebraic Analysis of Combiner Structure 





Let us consider what locking conditions, and hence what Boolean functions, are realized by the simplest 
combinations. 


A rigid element Y of any shape is controlled by interlocks from elements X,, Xzg,...X;p, and is free only 
for one combination of them. Hence, except for unused combinations 


























T (a1 Xe, Zq,--- 2m) = G(X, Ze, 
For oo p Ben) SPM Ze Meq.ss Lm; 
Ae ee | 
m4 y where xy = x; is set in this unique combination, and cro if xj remains 

k in its normal state. 

nko in the rigid kinematic scheme no displacement can occur if even 

4% % Im one element is locked. 
Fig. 5 


The situation is different if the locked element contains elastic 
couplings. Such a system is that of Fig. 5. Here Y has several (say, k) 
modes of motion, and so Y, controlled by locks from elements 
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X,, Xg,...X__ is free not just for one combination, but for k such. Hence the locking condition is expressed 
by the function 


} (X15 Le, Lg,--- Lm) = B(2y, GLa, ~~. 52m) = 
= (2:2 :9%13 see 1m) + (Z_1Ze0 see Tam) + (1) 
ES k kom 
+ (LeiFeg. ++ Tim) = Be (LisLig Lig ... Lim) = > I] (X43) - 
i—j] i—1 jel 


In the general case the combiner may contain several stages, as shown in Fig. 6. 


If stage 1 realizes functions fyy, fxg, fyg,---sfy),. 
the second fay, fze,.-+fgk, and the last (r th) stage 
fr, frgr+--efrk,» the whole system clearly realizes the 
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Fig. 6 If k = 2 the hinge system of Fig. 7 can be used 


instead of elastic couplings. 


Other ways of combining the elements can be analyzed similarly, In general, rigid couplings corres- 
pond to Boolean multiplication (logical operation ‘and"), and elastic or other nonrigid couplings to Boolean 


addition (logical operation ‘or"). Inversion (logical operation ‘no") corresponds to locking for other states of 
inputs, 


If there are several outputs, each has its own 
Boolean function to represent its locking condition, 


: : : 4 | In this way the structures, and hence the func- 
, tioning, of very complex combiners can be analyzed, 
I }— b The locking conditions for all outputs are written as 
Boolean functions appropriate to the coupling involved. 
To determine the output states from the input ones it 
Fig. 7 suffices to calculate the Boolean functions appropri - 


ately. This rapidly and exactly determines which of 
the outputs are locked and which free in any input state. 









































The formulae can be simplified, and hence the structures of the mechanisms improved, by using alge- 
braic transforms. 


Example. Consider the mechanism of Fig. 8. We write down the condition for element X controlled by 
interlocks to levers A, B, C, D, E and F to operate; 


F x = a6ee +020 + ase + 00 + a0e- (2) 


The transforms possible using Boolean algebra give 


Fy = 6 (a6e+ 20 + ce + 9) + 20e = 6 (ae + 9) + ede. (3) 
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The structure of the mechanism constructed from (3) is shown in Fig. 9. 
Different structures corresponding to these functions can be derived, and their equivalence established. 


Thus, for instance, if we transform (2) by a different regrouping, we get another structure 
Px =0(6,+ 6 + 0¢)+-a0e (6 + 1) = A(e-+ 2¢) + ase, (4) 


which satisfies the same conditions. 


Figure 10 shows the structure deriving from (4). Hence the mechanisms of Figs. 8-10 are precisely equi - 
valent, but the two latter are simpler. 


Algebraic Synthesis of Combiner Structures 





Use of the relations between ‘and’, ‘or’ and ‘no’ and the Boolean operations of multiplication, addition 
and inversion gives us the locking conditions for the outputs in algebraic form; these are usually given verbally. 
The Boolean functions are transformed to a form convenient for realization, e.g., that of (1), using identical 
algebraic transforms.* The general system of the mechanism is then easily sketched from the correspondence 
between the structures of the functions and the type of coupling between elements. 


It is sometimes convenient to use the somewhat artificial device of putting the operating conditions as 
an algebraic expression, multiplying the Boolean functions for the outputs by the symbols for the elements and 
adding, as in the theory of relay contact circuits ([ 4], § 1.6). An example is given below. The plus signs be- 
tween the Boolean functions are of logical and not structural significance, i.e., must not be interpreted as im- 
plying elastic couplings. 

Any Boolean function can be represented in the disjunctive normal form as of (1) [5]. Hence this method 
is always sure to produce at least one solution. The number of locks is minimized by minimizing the number 
of letters in the Boolean functions, and the number of elastic links by minimizing the number of plus signs 
( Boolean additions). 


Example, Let us consider the decoders used in telegraphy. 


We use the artificial device indicated above and write the operation conditions as a single algebraic 
equation; 


F = labcde + 2abcde + Sabcde + dabcde + Sabcde + Gabcde + Tabcde + Sabcde + Sabcde + 10abcde + 
+ llabcde + 12abcde + 13abcde + 14abcde + 15abcde + 1éabede + 17abcde + 18abcde + 19abcde + 
+ 20abcde + 2iabcde + 22abcde + 2Wabcde + 24abcde + 25abcde + 2Wabcde + 2abcde + 2abcde + 

+ 29abcde + 30abcde + 3iabcde + 32abcde 





Here the figures denote the letter levers. The letters denote the corresponding selection rods, those with 
signs over denoting positions allowing motion, the others preventing it. 


*The various identities of Boolean algebra are dealt with in § 1.7 of [4], and also in [5]. 
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In existing equipments all these combinations 
are distributed over the selection rods in accordance 
with the code used and with the positioning of the letter 


























| , | 
° levers. We neglect the distribution actually used, and 
fame Iq x transform the expression analytically. We get 
ie F = a(b] c [d( le + 2) + d(3e + 4€)}+ c{d( 5e + 6€)+ 
ey we | + d( Te + 8)} }+ b {cl d(9e + 10€)+ T+ (11e + 12e))+ 
+ c[d(13e + 14e)+ d(15e + 16e))} + a(b {cf d(17e + 
+ 18€)+ d(19e + 20€)]+ c[d( 2le + 22) + A 23e + 


+ 2de)}} + B {cf d(25e + 268) + U27e + 286))+ 
+ G[d( 29e + 30e)+ A Ble + 3ze)) } ). 


Figure 11 shows one possible design for a decoder 
to realize this expression. 


ys In this system the rods rotate and do not move 
axially. 


The normal position allows motion, the set 
position prevents it. 


In conclusion we hope that this exposition of 
the basic theory of combiners using Boolean algebra 
will stimulate further developments, and the intro- 
duction of scientific methods into combiner design. 





A The author is indebted to G. N. Povarov for 
12345 3132 valuable discussions. 
Fig. 11 
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THE TRANSIENTS IN A MAGNETIC AMPLIFIER CONNECTED VIA 
A RECTIFIER TO AN INDUCTIVE LOAD 


L. V. Safris 


( Rostov -on-Don) 


A method of calculating the duration of transients in a mag- 
netic amplifier connected via a rectifier to an inductive load is con- 
sidered. 


INTRODUCTION 


The transients in a magnetic amplifier connected via a rectifier to an inductive load (as in Fig. 1) have 
already been considered [1]. The deductions made there [formula ( 19)] show that the transient duration is 
determined by the load time-constant, no matter what the signal. The magnetization and rectifier charac- 
teristics were assumed ideal, with WL, >> Ry (where Lj andR; are the inductance and resistance of the load), 
and that the control current changes stepwise. 


It is of interest to see how far this deduction can be applied if the inductors and rectifiers have nonideal 
characteristics, 


Experiment indicates as follows. 


At control currents Ig close to maximal the rate of rise of the load current i; is determined quite accurate- 
ly by the load time-constant r; = Ly (Ry + 2r,) where r, is the forward resistance of the bridge rectifier. This is 
confirmed by Fig. 2, which shows an oscillogram taken on a test model (A, appendix 1) with 7 = 5.8x 10% 
sec. 


When the control current is reduced the transient is briefer. This is confirmed by Fig. 3, which is for the 
same amplifier at a control current 10 times lower. Curve 1 of Fig. 4 (amplifier A) was drawn up from this 
oscillogram. For comparison, Curve 2 corresponds to 
11,= 5.8 x 10” sec, 


sa The considerable divergence is explained by the 
1 back-emf induced in L; by the current change, e,; this 
h Wy acts to increase the voltage ujy acting on the rectifier 
iN, iy bridge. This can occur by the input voltage Uj being 

redistributed between bridge and output winding if the 


f R 
gt Way a P amplifier is not completely saturated (E,,# 0, Fig. 1). 
fe KA 




















ri Increasing uj,y reduces the transient duration (some as- 

tw pects of this have been dealt with in [2], p. 669). The 
‘saturation angle‘ will not be constant throughout the 
Fig. 1 transient, even if the control current changes stepwise. 
This is confirmed by Fig. 5, which shows u,,. The in- 
ductive load also influences the control circuit transient. 
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Here we deduce approximate expressions to be used in transient duration calculations for the inductors 
and rectifiers which have nonideal characteristics, and which allow for the control circuit transient, We use 
the linearized magnetization characteristics method [3], Only ‘weak signals in a choke amplifier are con- 
sidered. By ‘weak’ are meant signals for which the current and voltage acting on the working windings in no 
case pass outside the rectilinear range on the magnetization characteristics, so the amplifier characteristics 


can be taken as constant, The bias magnetization is effected with an extra winding (not shown in Fig, 1) in 
series with adequately high resistance. 
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1. The Equivalent Circuit 





It has been shown [3] that the ‘instantaneous mean‘ a-c current change in a resistive load Ai,,* (i.e., 
the arithmetic mean over one cycle recorded by a detector) for a choke amplifier in its linear range is 
RS 


Ai, = R,+R Ai, = k,Ai,. (1) 


Here kj is the current gain, Ai, the control current change (‘instantaneous mean’ neglecting even har- 
monics), R is the total resistance (sum of the load resistance Ry and a-c winding resistance rj; 1; << Ry, nor- 
mally, so in future we assume R = R,), S is the amplification factor, and Rj the internal resistance of the am- 
plifier. 


The amplification factor is given by 


where W, is the number of turns on the control winding ( Fig. 1), W,, is the number of turns on the a-c wind- 
ing, o = 0H,y,/OH¢ at B,, = const, H,, is the ‘mean‘ value of the a-c magnetic field component, determined 
from the ‘mean‘ value of the current: Hy, = W,1.,,/1 , Hg is the constant magnetic field component and is 
the length of the magnetic circuit. 





Fig. 5 


The internal resistance is 


Here Q is the cross-section of the two cores, yj = 4B,, / H,, at Hg = const, B,, is the ‘mean' value of the 
a-c component of the induction, determined from the 'mean' value of the emf applied to the control wind- 


ing Wy: By = Ew /wW_Q. 
uy and o are determined from H,, = f,(Ho, By), S and Rj as above, or directly from I, = fx{ Ig, E.,). 


o and yj define the magnetic properties of the inductors for a given shape of induction curve. If the 
range of variation in Ho, Hy, and By is restricted they can be considered constant. Hence S and Rj are constant. 
Thus for Fig. 6 (amplifier A) we can take S = 1.05, Rj = 3500 ohms (dashed curves). 


* The A sign indicates that the quantity can vary over the rectilinear range, and can be determined as the in- 
crease over the initial value. The mean i,, during the transient is a time-function. This is the sense assigned 
to instantaneous mean, as contrasted with the steady-state Al,,, 








It has also [3] been shown that, no matter what the load, the control winding inductance is 








= Ws Wav @(Ae,,) 
Hoa rete Ragen  — We d (Ain) * (2) 


Here c= 0By /dHg at By, = const, v = — OBg /AB,, at Hg = const, Le is the control winding inductance 
at zero load resistance and L,, is the additional control winding inductance due to the load reaction. 


Ht. and v are determined from Bg = fx( Hg, B.) 
and are taken as constant. 


The equivalent circuit shown in Fig. 7 has been 
drawn up from (1) and (2). It is here supposed that 
the delay introduced by the a-c circuit can be neg- 
lected and hence the emf is in phase with Aig, This 
assumption is based on the fact that for a resistive 


aa 
load the ‘instantaneous mean’ At, = RoR is al- 


most in step with the instantaneous control current Aig. 


In deducing (1) we assumed that the bias produced 
by the control signal Aly reduces AE, by an amount 
equal to the increase in load voltage, RAI,,. 


The experiments show that this remains correct 
in the steady state for an inductive load connected 
via a bridge if we take R as the input resistance of the 


Fig. 6 bridge Ry = 











UM 
Aly If E., and Ujy were to add arith- 


metically, as in ideal amplifiers [1, 4], they would be 

related by the straight line AB of Fig. 8, and so AE,, = 
=— AUjyy. The curves representing the relation of AE, to 4U}y for real amplifiers with inductive loads de- 
viate somewhat from AB, but have the same slope over a fairly wide range. Curve 1 of Fig. 8 (amplifier A) 
shows E., = {¢ Uy), as an example. It corresponds to the load characteristic of Fig. 6, which can also be 
approximated by a straight line. 





We here assume that in the steady state 


AE,=—AUy= — Ry Ag. (3) 


R Since the magnetization curves scarcely depend 
on the load type, the equivalent circuit of a choke 
magnetic amplifier with an inductive load can be de- 
rived from Fig. 7 by replacing R; by Ry, at least as 











regards the steady state (Fig. 9a), and so AI, = 
ee 

Ri + Ry 
replaced by a generator of emf RjSAlg and internal resistance Rj loaded by a bridge rectifier. 


Al,. In other words, the output circuit is 


Figure 9a has been demonstrated correct for the a-c circuit for the steady state. For ideal amplifiers it 
has been shown [1, 4] that the d-c and a-c ampere-turns are equal whether the load consists of resistance with 
or without inductance. This follows directly from the expression for Al, ,, if we remember that ideal charac- 
teristics correspond to 44 =@, o=1. 
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We now determine whether this equivalent circuit 
can be used to determine the loac current in the trans- 
{ent and steady states, provided the signals are weak. 


To this end, using Fig. 9a, we determine the cur- 
rent Ai; = f(t) and compare calculation with experi- 
ment, The transforms we shall use are correct if the 
final results are confirmed by experiment. 


We first simplify the circuit using [5], according to 
which the mean current in a (Ly,,R,) load connected 
via a rectifier to an a-c generator of interna] resist- 
ance Rj subject to a stepwise change in emf of 2E,/# 
can be determined from 

t 
gi ee | (1 #e), 


t RR, +R, (4) 


where 


Ly 
C=. 


The factors —€ and p allow for the effects of com- 
mutation on rectification, Figure 10 gives their values. 
these were derived from graphical calculations for sets 
of particular circuit parameter values, The effects of 
the forward and reverse resistances r, and r,, of each 
arm are allowed for the recalculating R; and Ry. 


The load current in our case is 


(5) 
=k, AI,(1— e *e) =Al,(i—e *). 


s' LJ 
Here k; = eee. is 


Ri + RE the current gain, 


= (R; + 7,) 0.5 , ’ : 
N= Ft” p= +h, S = 





R; , L 
Sa, = 
Ly 
oo” 
ec RAR, 


is the equivalent time-constant of the output circuit, 








Eq. (5)shows that the transients in the output are just those in a circuit consisting of an inductance Li 
and resistance Ry +R} on connecting to a d-c source, 


If Aly is a known time-function, our assumption that S, Rj, #¢ and v are constants means that Ai, can 
be found by Laplace transform methods, 


We can therefore replace this type of magnetic amplifier by the equivalent circuit of Fig. 9b. Of the 
two outputs (a-c and d-c in Fig. 9a), only the one to the load is here retained, as this is of direct interest to 


us, Hence Fig, 9b is equivalent only as regards the load circuit, and gives no indication of what occurs in the 
a-c circuit, 


To find the accuracy and applicability of (5) when 
deriving the steady-state Aiy, we use Fig. 11, which 
shows 1) I; = fx 19) (amplifier A), and 2) the same cal- 
culated from the graph proposed in [4] (Fig. 4) for an 
ideal amplifier, plus 3), which shows Aly, calculated 
from (5). In the range Ig = 2-20 ma (5-50% of Ignax), 
which is of greatest interest to us, all three curves prac- 
tically coincide, Near the maximum I, the divergence 
is large, because we have taken € as constant, The 
measured and calculated Al, for other cases are given 
in Appendix 2, 


We now consider the transient duration, 


2. The Transient Duration 





Three particular cases are considered, 





5 0 tH 0 HB MM & 


Ai, varies stepwise. The tansient duration will be 
determined by the time-constant of the output circuit 
Te» 1.e., by (5). From the latter we have drawn (Appen- 
dix 2) the Aj; /Al, = fg(t) curves for the models, A, B 
and C which differ in the inductors and rectifiers. These curves are shown in Figs. 4, 12 and 13 (curve 3), 
Curves 1. on these Figures were drawn up from Figs. 3, 14 and 15, which were taken under conditions corres- 
ponding to those used in the calculations, The divergence between curves 1 and 3 (at Ai, /Al, = 0.63) is on 
average 30%, and arises partly from the finite rates of change of Aig and ABg, and partly from errors in mea- 
suring the oscillograms. In one case (Fig. 4) the rectifier back resistance is neglected, 





Fig. 11 
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Fig. 12 
The time -constants given by the oscillograms are several times smaller than r (curves 2 on the figures), 
and differ more from one another the larger R} /R; 


Aug varies stepwise, Here the law followed by Ai, is determined by the parameters of the output and 
control circuits. 











The equivalent circuit of Fig, 9b shows that the 
transient on the output is described by 


d(Ai;) (6) 


S Riad, = (Ri + Ry) Ai, +L, —T EY 














The equation for the control circuit (using (2) and 

















Fig. 9b) is 
. . W2Q acai Wov @(Ae,,) 
Au, = R, Ai, + Uc ee 7 Wo at « (7) 






































The Ae, in the ac coils during the transient we 
take to a first approximation to be 








d(Aiy) (8) 


— Ae, = Au, = (Ry + ra) Aipt+ Ly dt 





Unlike in (3), we here assume that RMAI,,, = (Rf + r,)Al;,, and also that the emf produced by the load 
self-inductance from Ai; is equivalent to an equal increase in the voltage on the rectifier terminals, This 
assumption is requirea because the output circuit in Fig. 9b cannot be used to determine ey, = L/d(Ai, ) /dt 


([6], p. 48). 


This simplification is justified by the results from the following experiment. When an opposing emf 
Ey, ( 0-8V) was connected in series with the load (amplifier A) E,, took the value given by curve 2 of Fig. 8. 
The E,, = f;(E,) calculated from (8) is shown by curve 3, The divergence {s slight. Hence using (8) we allow 
of a certain inaccuracy in determining the effect of the load on the control circuit transients. 


Applying the Laplace transform to (6), (7) and (8) solving we get 

















. AU 1 i 
A oak coe - —1 a 
rh. Ro kL 4 FP TEEPE |’ (9) 
Ai, = Ser (= Tie 
Pa P (et) + ++ tpt! (10) 
Here 
: _ Wo Qe a Wov (Ri,+ ra) ky ae Ly (11) 
¥ SE W ,, oh Jaa ° Ritts 


These expressions indicate that the control current does not rise exponentially during transients, This is 
confirmed by the Aig oscillograms of Figs. 16 and 17, for amplifiers A and B, with Aug changing stepwise 
( Appendix 2). 


Applying reverse Laplace transform we have 





_— (ett WAV te +t, FTP — Ft + V5) (12) 
aad gen 2(ToTe+ TT) 














Fig. 14 





























Fig. 16 


Two main cases occur; p; and ps are real (—a,,—ag) or complex (—a2jb). 


In the first case Aiy is 
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Ai, = 











Fig. 17 


If the roots are complex 


: AUg 1’ b 
Ait = 2 ki [1 — Fe" sin (Ot —g)], - 


where 
1 


ab wan ay EE ee 
9 = arcan— , B= a+b? = ert, 


Hence with an inductive load connected via a rectifier we can get damped low-frequency oscillations 
on Airy, 


(Fourteen) has been used to determine the Ai; /Al,= f,(t) shown in Figs. 4, 12 and 13 (curves 5) for the 
A, B and C amplifiers (Appendix 2), They correspond to the oscillograms of Figs. 16, 17 and 18 from which 
the curves 4 in Figs. 4, 12 and 13 have been construc- 
ted. In all cases the curves 4 and 5 coincide adequately. 
The oscillatory component appears more clearly on the 
calculated curves than on the oscillograms, 


The control current varies exponentially. This 
case is of some interest in relation to approximate 
estimates of the lag in the current in an inductive load 
from the lag in the control circuit, if known. 





We assume the control circuit parameters R, and 
Lg constant, Hence Al, in response to a step voltage 
input will rise exponentially. 


The expression for Ai; can be derived from (9), 
(12) and (13) by replacing r . by rg = Lg /Rg andr, = 0. 


Substituting in this way, we have 





, AU, a ks 4 
dip= zthl— pasa Haren): 

















The curves 6 corresponding to this case are shown in Figs. 4, 12 and 13 (Appendix 2), There is no great 


divergence between the results given by (14) and (16), We must remember here that r 4 was taken direct from 
the At, oscillogram as the time corresponding to Vie /Aig = 0.63, 


SUMMARY 
1, The equivalent circuits and formulae used are confirmed by experiment. 


2. The expressions for Aly, indicate how the amplifier and rectifier parameters affect the transient dura- 
tion for an inductive load and weak signals. Estimates from the load time-constant alone are useless. 


3. The expressions are suitable for use until stricter and more universal methods for the purpose become 
available. 


Appendix 1 





Amplifier Descriptions 











Model 
Data A | B C 
Cores Material EN80 MZ Permalloy Transformer 
(NIITVCh) steel 
Construction Spiral toroids Spiral toroids | Ring toroids 
10-*Q, 20.325 2x1.5 21.67 
1, 0.13 0.2 0.25 
Windings W..» Turns 190 40 800 
Ww, » 200 120 600 
r,, ohms 2x2.4 2x0.8 2x 2.05 
Parameters Ss 4.05 3 0.55 
R,, ohms 3500 105 1600 
10-* ,, 0.6 3.8 2.5 
v 1,38 1.35 1.5 
Bridge Type DGTs-4 Curpaus oxide | Selenium 
reotifier Tr,» ohms 122 9 40 
» » co 2200 9000 
320 
Ly, Henries 25 19 . 
Frequency, c/s 500 500 
ae U,,¥ 16.2 22.5 63 
conditions L’ 
lop (bias), ma 5 15 45 




















Appendix 2 


Calculation of At, 

















Model 
Operating conditions R " " 
R,, ohms 3620 103 1200 
Ry, , ohms 308 109 360 
Ri Joly, 0.004 0.0015 0.06 
9 ls BF 8 pe oe ie 
ee when! : : : 
Sig= | Alowhent>0 pe 0.63 0.93 0.8 
Al,, “a 5.0 30 30 
exp. 1.2 13.5 1.9 
tor sec. | OTe. 0.8 (67%) 9.6 (71%) 1.5 (79%) 
Ai; | &*P- Fig. 4(curve 1) Fig. 12( curve 1} Fig. 13 (curve 1) 
ZT, } caic Fig. 3. Fig. 4| Fig. 14.Fig.12 Fig. 15. Fig. 13 
L . (curve 3) (curve 3) 
(curve 3) 
Ry, ohms 6.8 4.5 134 
10-*t,, sec. 0.18 1.8 0.9 
10-7, > 3.65 4.85 0.52 
A { 0 whent<0 10-*¢ sec. 5.8 16.1 4.75 
“o= | AUWwhem>0 po, 1/sec. |—10.84/18.5 | —8.54/5.7 | —38.3/ 34 
tan? —1.71 —0.67 —0 885 
?, rad —1.04 —0.59 —0.725 
by1/sec, 21.6 10.3 51 
bo/b 1.17 1.8 1.5 
exp. 7 46 14 
Al, .me calc, 6.3 43 14.5 
Ai exp. Fig.4 (curve 4)| Fig.12(curve Fig. 13 (curve 
are ss Fig. 16 Fig. 4| 4), Fig. 17. 4), Fig. 18. 
c. 5 ig. 12 Fig. 13 
. mea! feive 5) (cirve 5) 
, 10-*t» sec. 7 7 1.6 
eg a 4 { 
ai- = 1—e * 402 0.16 0.38 10 
e Te—-To’ sec. 
Ai, ; Fig. 4 Fig. 12 Fig. 13 
va calc. (curve 6) (curve 6) (curve 6) 
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THE CALCULATION OF A CHOKE MAGNETIC AMPLIFIER CONNECTED 
TO A SINGLE-PHASE BRIDGE RECTIFIER 


N. A. Kaluzhnikov 
(Kharkov) 


A method far the calculation of a choke magnetic ampli- 
fier connected to a single-phase bridge rectifier as regards the 
a-c circuit is presented, for the case where the load has induc- 
tive and capacitative components, using the theory of the ideal 
‘linear’ magnetic amplifier. 


Practically the only papers on the calculation of a choke magnetic amplifier connected to a single- 
phase bridge rectifier are those of Komar [1] and Vasil'eva [2]. 


The calculation of a complex nonlinear circuit such as the ac circuit in a magnetic amplifier connect- 
ed to a single-phase bridge rectifier must involve some fairly major simplifications, 


Here we also make some assumptions, 

1. The rectifiers are ideal (forward resistance and back conductances both zero), 

2. The ‘natural’ magnetization mode is used (impedance for even harmonics zero). 

3. Core losses are neglected, and the resistance and leakage of the windings are assumed zero, 


4. Only large inductances and capacitances are considered so that the load time~-constant is much longer 
than one half-period of the supply. In other words we assume the inductance to smooth the current ideally 
and the capacitance the voltage, These conditions are adequately satisfied in practice by a load with w1) = 
(5—10)Ry and wC)=(5— 10)/R. 


Here we consider two extreme cases, 
1, The cores consist of a material with a rectangular magnetization curve. 
2. The magnetic amplifier acts as a linearly varying inductance (linear case), 


The latter case is only of interest in calculations as representing an extreme case, This type of assump- 
tion is made if the currents and voltages are assumed sinusoidal with the load operating on a-c, which is the 
basis of most existing calculation methods, 


Real amplifiers work in modes lying between these extremes, We therefore need only consider these 
modes, and recommend the less favorable of them, so that normal control operation is possible under real 
conditions. 


All load parameters (resistance, opposing emf etc.,) are assumed given, together with the bridge output 
current I). 


The object of calculations on the a-c circuit is to determine the supply voltage U,, the minimum voltage 
on the amplifier Ua, and the current drawn by the a-c windings Ig. These are all calculated from the load 
parameters and from the magnetic conditions used, which are determined by kg equal to the ratlo Us max /Ua, 
Us max is taken as equal to Uc. The calculations are done using half-period mean voltages and currents. 
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All types of load are divided into two groups, namely Komar's A [1], where the impedance {s inductive 
and ‘double short circuit’ operation occurs and R, where the impedance is capacitative and ‘double open cir- 
cult" operation occurs, 


The purely resistive load case falls between these and allows neither A nor R modes. 


Operation with an Inductive Load 





Figure 1 shows three circuits (a, b and c) to be considered in this section. They are ones with an emf 
acting in the same direction and a resistance, with an inductance plus resistance, and with an inductance plus 
two resistances, In all cases there is a divergence between the mean currents before and after rectification, 























Fig. 1. Circuits using single-phase bridge rectifiers and inductive loads, 


Ideal amplifier. We start with the circuit of Fig, 1a. Figure 2 shows graphs of ug, Uj, la and 1). 


The operation of the circuit falls into two sections. In the first no current at all is passed, and the line 
voltage is entirely taken up in magnetizing the cores. During this period a current 1; = EgR)flowsin the load 
circuit, 


During the second the amplifier is completely saturated, 


The circuit then presents zero impedance, and the rectifier begins to conduct at the instant of transition. 
The total current flowing is 1; = (uj + E9)/R;, so uy = 1) Ry — Ep, 


A more detailed description of the operation is given in [1, 2]. 





The mean values of these quantities are 


U, =+ Um C8 cot deot = — tm (sin) + 1) =U, +S | (1) 


—n/2 





1 











n/2 
i 1 : i—sinyg 
U; = = Um cos at deat = = Um (1 — sing) =U, —,—_, (2) 
{ i E {4 u i ¥\E 
ras, 2 m fo as _m —_ gi |e 
igs \ (a coset +p.) det = m Ry (4 -sin #) + (5 DF tes 
Y 
4 U. i—sed (5 ) '\ Bo 
= Ry 5} - I~ x)R,’ (3) 
d n/2 : 
eT Ey 1 ( (Mm Be, , Ue 1 ee Be 4 
te | dat + \ (qr cost + Ry) dat = ainhe b ” 
—n/2 y 
The expression for the two coefficients used are 
: v, . (5) 
P Me i+sing 
and U 
(1 — sin y) { y ‘4 Y\ Zo 
c 2. 1+ r->)T 
pate eee 
i—singd ‘ns E. : 
1 o—5— + Be i+7, (6) 


Here 


Uy = hy —E, =U, +See 


We turn now to Fig. 1c, which is the equivalent circuit with allowance for the mutual inductance (Fig. 
1d) at unit coupling coefficient (M*/LiL, = 1). The parameters of the two circuits are [1] related by; = R;, 
t= Rg (Ly/Lq), 1) = Ly. 


Using Komar's results, we have #1=1,;R—Ijr where R=R, +1. We thus get a circuit equivalent to 
that with a forward emf, if we put Eg#I,r. This enables us to transfer the deductions to this circuit also, The 
two circuits show an essential difference in mode of operation, however. In Fig. 1a the emf does not depend 
on the magnetic amplifier at all, whereas in the latter case the emf depends on the amplifier current, 


Hence Ug, Uj and kg are given by (1), (2) and (5), The expressions for Ig, lj and8 are somewhat altered. 


Since the d-c component of the voltage across the inductance is zero, we have 


For I, we then have 











Substituting for I from (4a), we get (after manipulation) 





= 
Po ace Semee (yf 4 oo) er 
—— * [1-( x) +r} (3a) 
+ 2) r i 
p=1—— = r— 2p re 
E 
ho ae 
(Six) and (6a) are identical if Ui = Rj = kp. (6) 


Hence the graphs derived for Fig. 1c can be applied to Fig. 1a if the kp are equal. Suppose now that kp 
varies from 0 to a due to change inr at Ry = const. 


For kp = 0 we have a purely resistive load, and (1), (2), (4a) and (5) apply. The expressions for I, and 8 
simplify to 


7. — ve 1—siny 
ents 


‘ =J,and p=1, 


For kp = @ we get the circuit of Fig. 1b. Here 





7. — Ye 1—siny x—2y 


a Rj 2 an (3b) 
_ z—2h 
p-=—- am 


The analysis therefore gives this result. For all inductive load types without exception U, + Uj = U, if 
the core material is ideal. The extreme cases are kp = 0 (pure resistance load) and kp = @ (resistance + in- 
ductance}. In the first case 6 = 1, in the second 6 is maximal. All other cases are intermediate. The non- 


linear distortions introduced by the load are defined by kp, which is r/R for Fig. 1c, and E9/U, for Fig. 1a. 
Figure 3a shows 6 = f(kg) for an ‘ideal’ amplifier. 
Linear choke. A bridge fed via a linear input inductance is treated in the same way. 


Figure 4 shows how the currents and voltages in the circuit of Fig. 1b change in time. Unlike in the pre- 
ceding case, the two halves of the circuit operate independently in the first section. Here 


: 4 E, : u u 
= Ue, Uy =i, R,— E,=0, iy= Zr andig = \ — dwt+C= — sinwt+C. 
oie Gabe | 0 i> R a \= + X, i 


di 
In the second section the law followed by the current is given by u, + Eg= X, ry + igR,, where i, = 1), 
di 
uy = i,R; — Eg and Ug = Up ~ Wy = Xg —, 


These give the instantaneous currents and voltages as 


. 
t 


u 
is = iy =F cos (at — 9) + 7 + Ae 
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xX 
(where Z = / R? + x , tan g=T= — , and Xq is the linear input inductance); 


u ot 
u, = 7 Ri cos (wt — 9) + ARe * 


* wt 


Xx os 
sin (wt —¢)—-—&-Ae *, 


um 


u =-—X, 7 


E E 
The integration constants are given by the initial conditions ig (@) = zh and i, (¥)= R . 


This gives the choke current variation law as: 














E 
in the interval (y— ©) to @,ig, = va (sin wt — sin @)+— ; 
u sf E 
in the interval @ to y, 1, = ij = _ cos( wt g) —cos(@ — g)e : + ;"? 
in the i Lee = oe (sin wt at a 
n the interval y to(@ + mig, = X. (sin wt— sin p~)+ . 
Expressions for u, and u, in the interval from @ to are represented: 
~ S53 
1 =,tm.008 9 [cos (tat — 9) — cos (9 — 9) ¢ * 1, 
~& 3 
u , = —um|sin gsin (at — 9) — cos pcos (6—g) e * |. 
The formulas for the mean currents and voltages are 
E if E E U. sin —sin® 
. c — 
hag ts\(i—#) eet ea. (7) 
8 
E if E 
. 0 
Is =H +5 \(i-H) dwt -} 
if E 1 “é*, E 
) . 0 
+ =| (i,, —F) det + — (i —F) dat = 
U. 1f0—a .. —*+.. cos — cos —2 
-- =| >= sind + =t4 Y sing + a =e), (8) 
aS _ yy sin? —sin 6 9 
1f i : i 
U, = =| we deat + = \u dot — 1 np dt — 
a 6 
oe 0 0 
—— \ ue dat = U.[+cosd-+ 2 FS2"), (10) 
Y 
ipa 3 11 
B 2¢cos 8 + sin + 8in®’ (11) 





























U ne ed . 
a +(° @ sing +2 +4 ¥ sin » + 289 ons ¥ Reo) 





- Ey 
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Ue sin » — sin 0 


-}- 
1. 4, 2 (12) 





E 
The angles », 5 and a are given by |; (~) = 4 
ug (5) = 0, lay (a) = 0. ; 























u The resulting transcendental equations for the 
. angles are : 
| 
+ ‘ o—0 
Zz Hed 1 ae IM, on Soa 
Z ly oe wt cos(}— 9) =cos(9—)e =, (13) 
2 
U 
1 
3-0 
a; / i wt (6— 9) THE gre @ » (14) 
i 
t 
\ \ sing == SUO—sing (15) 
Lt ties , 
? Es ay dn 2 wi In deducing the last equation we have also used 
E E 
‘a the condition i, (@ + ™)=— = , which gives — = 
1 l 
x i u 
te i = = eas . This latter, when transformed, 
$< wt 2 Xa 


gives a relation important to the subsequent analysis; 











Fig. 2 oie ae, HOY ohm 
af 2 Ey 2 ered 
Substituting for Uc from (9), we finally get 
ae U, sin }+ sin6 
2 Fo/ sind — sind” (16) 


Equations (9) — (11) and (13) — (15) continue to apply to the circuit of Fig. 1c. 


(Seven), (8), (12) and (16) are somewhat altered: 











=“ Ui U. sin ) — sin 6) 
af “kane al 2 ): ( Ta) 
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Fig. 3 
Rh 1;0—a r+a—y cos 8 — cos b — 2 cosa 
we 1 r+a—y . 
Ry+r tl 2 in + 3 o + 2 ] 
P=i>= sin} — sin® ’ (12a) 
2 
gene A ein) + sind (16a) 


2 r sing’—sin0 





X 
Also rt = R om . This gives the required relation between the circuit parameters and @ and ». 


Let us consider the limiting cases r = 0 and r= @. 


When r = 0, we get on passing to the limit that 


U. 
U, =U, cose, Us =Ucsing, 1i=/a = p-cospand B= 1. 


We consider the second case, which corresponds to the circuit of Fig. 1b, in somewhat more detail ( Fig. 
5 illustrates the processes occurring). 


Here r = Oand g = 0 whenr =. Now (13), in which the RHS becomes zero, gives cos(~ — y)= 0 
T 


agus 9. . Hence sin » = cos g= 1. 
( Sixteen) then becomes 


Xa = i+sin®@ 
She Ry, 2 i—siad (16b) 
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(r, is the time constant for the circuit of Fig. 1b). 


The other equations also become 














ae U. 1—sin® 

aes OR) Sa (7b) 
i —sin® 
| ek eee acre ( 9b) 
Va, =U- + (10b) 
ke = 7 sind llb 
+ sin 6 (11b) 
; 9 Oa x $a Ola It wt 





y 
r eS 1—sin® 
On iis a cows ei: (15b) 
«le aaa 


Fig. 4 





U. 1 1—sin0 
{ =hy— Fe TH sind [OO — 28nd + (5 +2) + cos 6 — 2cosa, (8b) 


(@—a)sin 0 + (> +2) + cos 0 — 20s a 
} i+ sin6 





p= 1— 


alr 


( 12b) 


Let us summarize the results. For a linear choke Ug + Uj = U, applies only if an inductive resistance is 
connected to the bridge output. The deviation from this law increases as kp falls, and becomes maximal for 
a pure resistance load, to which applies the iaw us, + Uj = us. 


As regards 6 , the position is the same as for an ideal amplifier, i.e., 6 falls as kp increases. 


The worst of the two extreme modes is, as regards the voltage equations, the ideal case, since Ug = U o Uy 
is the minimum U, which can occur in the circuit for given U, and Uj. In no case can U, be lower than this 
calculation gives in a real amplifier. A real amplifier will therefore not have a performance below the calcula- 
ted one, 


The error in U, is minimal for a load consisting of a resistive inductance (in both cases the law is 
U, = U,— Uy) and maximal for a purely resistive load. 


The 6 for a linear inductance was derived above. This has also to be estimated in the calculation, since 
for a given 1; the I, will increase with (1, = 61;). The most reliable calculation is that which gives the highest 
current for a given Ij. 


Figure 3b shows B = f(kp) for kp = var and a linear inductance. 


U 
Figure 6 shows the load curves in relative units for kp = 0, 1 and o, where lj ,,,, = rs 


Operation Into a Capacitative Bridge 





This type of circuit shows a divergence between the mean voltages before and after the bridge. Three 
such circuits, to be considered here, are shown in Fig. 7. 
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Fig. 6. Load curves for an amplifier working into a 


Fig. 5. bridge rectifier with an inductive load. 


As before, we start with an ideal amplifier. 























Fig. 7. Magnetic amplifiers used in operation into a capa- 
citative bridge. 


Ideal Magnetic Amplifier 





Let us first consider the circuit with the opposing emf and pure resistance (Fig. 7a). 


Figure 8a shows how ug, uy and i, vary. 


The operation splits up into three sections. In the first (from — 1/2 to — @) the rectifier does not con- 
duct because of the battery ( ‘double open circuit"). 


From—@ to the line voltage is distributed between the amplifier and the rectifier; no current flows 
because the amplifier is unsaturated (core flux increasing), so ug = Ug ~ Eg. 


uc —E 
From y to @ the amplifier is saturated and all excess voltage (u, ~ Eg) acts on Ry (1, = —" 


Uy = Up + Eg = iyRj + Eg). 








The mean values of I, U) and U, are 








6 8 
4 ( U4,» cos at — EL, _ 4 04m a 
I, =1, == R, fot = |g nse een fot 
oa 
U. ¢sing—siny O6—4% (17) 
= 7 ( 5 ——- cos 0), 
where 
obs 0 w-He. 2 & . 
Um n U, 


U; = 1,R, + Ey =1,R\ += U-cos8 = 








~ a oot (18) 
=U, (Se me ys g+* cos 0); 
‘ y 
U, = — | tm (cos eat — cos 6) dest = 
—@ 
= U,( ee rane 2+ 8 008 6), (19) 


The expressions for kg and k,, and one for Kp as a function of @ and », which is required later on, are 














oe 2 

B™ Ua ~ sing + sind —(o + 0) C086 ’ ( 20) 
<a : 

ce Uy ~ sing —sio > + (x —@ + 9) cos 6’ 

_ 2 6 2 £E, fsing—siny 6—¥ (21) 

cos 6 = heer if | 5 ig cos 0], 
so 
hea TR, _. sind—sind _6—y 
Eo ™ cos § cepa sy (22) 


Consider now the circuit of Fig. 7c. In this case [1] uy = ijr + 1yR. Hence all deductions for the previous 
circuit can be transferred to this one if E, is replaced by 1j)R and R, byr in (17) —( 22). 


Putting r + R= R, for Fig. 7c, we have 





es U, Ry /siné—sing 6—v 
h= 3 | 5) _ 5) cos 0), (17a) 


Since U, * }; (r+R) = 1yR, [1], then 


R ; ol 
UV, =U,—h(™mersing —* 





+ cos8), 
(18a) 








Equation (19) and (20) remain unchanged, while (21) and (22) become 








PE moeY 2 
1 y sin @ — sind —(6— $) cos6 ' “aia 
k _ Ji" _ + _ sind—sing 6—y 
: — 6, a cae (22a) 


The processes in the two circuits are the same if the kp are equal. In addition to this important con- 
clusion, which enables us to calculate both in the same way, we get a futher useful expressiont 


rst _ sind—sing 6—y 
7 00s 6 = z + cos 4. 





With this we can simplify (17a), (18a) and ( 21a); 





R 
a ee wel coe 
= + cos 6, (17a ) 
R nr 
U, = U.—} = cos 8, ( 18a") 
reat set ae 
"| o_o ( 22a") 
When r = 0 (kp = 0) we get the circuit of Fig. 7b, for which the formulae are 
U 
= m + cos 6, (1b) 
1 
U, =U. cos, (18b) 
U, =U, (sin 6 — 6cos 6) ( 19b) 
(the last derives from (19) by putting ¢= 6), 
i 
ke = Sa0— Oct’ (206) 
es er 
oe cond cae 
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Fig. 8b 





Figure 8b shows how ug varies; uj equals E, 
throughout here, Then i, takes the form of puises 
of infinite amplitude and infinitesimal duration. 
It flows at wt=@, +6, 2" + @ and etc. The 
mean value over a half-period is U, /R}. 





























When R = 0 (kp = 0) we get the pure resist- 
ance case. Hence Fig. 7b has the maximum non- 
linearity, which becomes less as kp increases. 


Fig. 9 


The result of this analysis is the k, = f( kp) 
curves of Fig. 9 for kp = var, with which the a-c 
circuit of a magnetic amplifier can be calculat- 
ed for given load parameters. 


Linear Inductance 





The inductance is connected to the input of a bridge which has a resistive load plus opposing emf, Two 


cases occur, namely,continuous and discontinuous current flow. The u,, u, and i, curves for both cases are 
shown in Fig. 10. 


The current is found by solving 


di 
a. re + Ri, =upcos wt — Ep. 


E E 
The initial condition isi, (—@)= 0, where @ is defined by cos 6 = — 5 ae . This gives 


i _ __ at+0 E _wt+6 
im = J [cos (wt —9)—cos(8+e)e * J—T*(1—e * ), 4 =I! 














P X 
Z= x? fo 2 dt = a 
where V & Ay and = 7 ae 


uy and u, take the forms 

















, ‘ an : , __ ot +0 at $0 
u,=i, R, + E,uyg =Xa oo =U» [— sin psin (wt — ~) + cospcos(6+ o)¢ * ]— E,e ? 
(for more detail see [5]). 
The mean values of lj, U, and U, are 
1 * U sin } + sin 6 y+ 6 
W=T, = = \ igdut = [ 2 2 cos 6], ( 23) 
—6 
| * in} + sin® ) 
U, = = \ u dwt = U. [= ~_ + T= $=" cos 0] ’ ( 24) 
—0 
, ’ 
Uz : = \ u gdut — = \ Ma deo = Ux (cos 8— cos 9). (25) 
—v & 


Here y and 6 are found from i, (~)= 0 and u, (5) = 0, which give us transcendental equations 























a 28. cos 8 ~2t* 
cos(p}—) =cos(6+)e * + apie FF), (38) 
8-40 
MS hoe, _ 40 
sin (8 — 9) = co Ot oe _ os “Rtn 
? (27) 
The expressions for kp: ke and kp are 
he = - 28 
Bb “= (cos 8—cos6) ’ (28) 
ak 2 
ec ~ sin} + sin ® + (x — »—8)cos6 ’ 
; Ft ee E. sin } + sin 6 _ +8, (29) 
cos § = — U. ast rk, ( 2 2 cos), 
so 
kp — Ji _ sind+sind  4+0 
 — — wee  . (30) 
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Fig. 10a Fig. 10b 











In the second case, when the current flows in the inductance throughout the half-period, the expressions 
change to 








; Um ~ S222 E _ at+n— 
iy = —z7- [cos (wt — 9) + cos (p— o)e . hw i (i—e an 
iy = ig |, uy = 7, R, + Zo, 

ot +n—y __ at+n—y 





Us, =Um[—singsin (wt —~) —cosecos(p—g)e * ]—Ee t 


These formulae ate obtained by putting 6 = 7 — yp: 


y 
1 Lt 
J, =], = — \ i, dut = 7~ (sin p — 5 cos8), (23") 
y-—r 
Uj = = ¢ ; ’ 
i= = \ u,dwt = U,sing, ( 24") 
y—r 


A ” 
U, =+ J ta dot —  \ugdat = Uc (cos — cost) 





$20 ( 25") 
The following formulae also alter somewhat 
a 
cos(p—g) = 20s te (26") 
cos @ a? 
ite * 
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so 











= Si n+3—¢ 
sin (6 — 9) = — say—9) , os ee ee (27°) 
tT sin @ ’ 
1 ( 29") 
ke = sing ’ 
_ 2sing 
‘n= xcos0 1 (30°) 
The boundary between the two cases is found’by putting » + @ = # (30), Then we have 
T “ 


Considering now Fig. 7c, and using the above expressions, we get the formulae applying to this case 


only: 
U. F {sind +sin0 y+ 6 
aga) 2} 2 gs cos 6], 











Ry [sin + si +0 
U; = hr+£,=Ir+->U,cos6= U. ! (and isne _¢ 5— cos), 














ke — 7. — Sinvt+sino _¥+8 (30a) 
a wane. — Elie 
The last of these gives 
si sin 0 - 
== mar cos8 = 1 ™ cos, 
so the expressions for I; and U, can be simplified to 
U. R 
meget. ies 5 
N= 7) HF 8? (23a) 
R 
U, =U. ~ cos 0, (24a) 
es: ae { 
he = ee eos (2%) 
Equations ( 25) — (28) retain their forms for this load type. 
If the switching is continuous the relationships also alter somewhat: 
U. R . U, U. R 
_ _1 [sj a 50) ax —a- n 6 en ot 23a 
I, = i, (sin }— -5 cosO) Ry sin» i, 2 cos 6, (23a") 
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U, =U, — =! (sing — = -;- COS 0) = 


cos 6 = 4 
ria 


ho = 


U.siny = U.— cos 0, (24a") 





Ry 


3 
2sin 
hae Bias. 
(30a *) 
i sin , 
1 oe 
cos6 sin} * ( 29a") 


We now consider Fig. 7b, whree kp = 0. Figure 11 shows how the voltages and currents vary in time; 


the relevant equation is 


di 
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Fig. 11 


U, and the coefficients are given by 
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= Um coswt — Ey, 


the solution of which, subject to the condition i, (-6)= 
= 0, will be 








. Um . . 
i, = x, [sin wt + sin 6 — 
— cos 6 (wt + 6)}, 
where 
J,R 
cos@ = —t-1 — _% 
Um Um 


The voltages are found from u, = u, —E» (u, = 
= Up COS wt — 1)R,) and uy = 1,R; = const. 

For this load type we get, without integrating, 
U) and I;. 


U,= 1, Ry =U. cos, ( 23b) 
U, 7 
I; =F 7 0086. ( 24b) 
6 
ae 4 
= — Um (cos wt — 


0 
— cos 6) dwt = 2U7, (sin 6 — 
— 6 cos 9), — 
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Rit { 
B® ~~ 2(sin0—Ocos0)* = 28b) 
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mcos@ ° 


ky on (29b) 





If the switching is continuous the expressions for 









































: d 7 f é W] U) and I, remain in force, while U, is found by inte- 
Fig. 12 grating: 
5 ¢ 
U, == \ Um (cos § — coswt) dwt = U, [sin 6 — sin b + (pb — 6) cos 6] ( 25b‘) 
6 
and 
kp = : ( 28b) 
sin 6 — sin } + (fb — 6) cos® * 
To determine y~ we use the condition that 
4 y 
uy Um 
—_ (cos wt — cos 9) dwt = = \ (cos 6 — cos wt) det. 
v7 8 
This gives sin py = . cos@ (> @). 
The boundary between the modes occurs at 
u a Fo 
_& > ~j 4— bed 
4, Sa "iy a oo tg — eX ° ( 31b) 
08 \ Sane 
N \ Yen’ Figure 12 shows the results as k,. = f( kg) for kp = var. 
\ 
“ NRT Sane Figure 13 shows the load curves for a magnetic 
\ xT X amplifier with a capacitative load, with kp = 0 (Fig. 
a = : nm 5 7b), kp = 0.1 (Fig. 7c) and kp = oo (pure resistance 
a BN s fy ‘ load), for an ideal amplifier and for a linear choke. 
tg= \ 
a =. Nv % i The curves of Fig. 13 show that the calculated 
Ba A | Ug is least for the ideal amplifier. These curves should 
0 therefore be used to give a reliable design, The extent 
a2 a4 Qe = 4 of the error can be estimated by comparing the curves 
1 max for a given value of kp. 
Fig. 13. Load curves for a magnetic amplifier These curves are drawn in relative units, 
with a capacitative load connected to a bridge uv, end _ where 
rectifier. The full lines are for an ideal ampli- U.. 1 max 
fier, and the dashed for a linear input choke, U. R 


U. r 
me (1 + kp) > (cos) max 
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The value of (cos ® )nax 18 given by (22) with » =—@ (i.e., the amplifier is assumed completely saturated 


throughout the conducting section of the half-cycle), This gives kp = = (tan 6 — @), which is used to cal- 


culate (cos 9 max’ 

— ee ey ery Ee 
. ure shows 

R, g Ww max Ry 





For a capacitative load (Fig. Tb) (cos 6) 
= f(kp). 


ul 
max = 14nd max = “> 


Sequence of the Calculation 





The following sequence of calculation derives 
from the above discussion. Ij, Eg and the load impe- 
u,/fy dance must be given. The first step is to determine 

\ Uj and kp. Then, using the magnetization curves (Fig. 
if 15),we select the two working points (1 and 2) which 
correspond to the maximum load current and to I, = 0. 








12 





/ (np) The first should be chosen at the knee in the 
Gy) ee 2 curve; the choice of the second (of coordinates Brin 
I 
105 F, 4 rT) 7; 01,1 H. max) depends on the purpose of the amplifier (see, 
7 for example, [ 4]). 
































Fig. 14 
Then U,, Ug and I, are calculated from 
U.= kU, (32) 
U, 
Adee “he (33) 
I, = Bl). (34) 
For an inductive load k, is found from 
ky 


( this formula is easily derived using 

U. 
U.=U,..+ U, = % U1) , 
k. is determined from Fig. 9 for a capacitative load. 


kp is calculated from 
kg = ee (36) 


where B_.,, and B_,,, correspond to the points in Fig. 15. 
For an inductive load we select 6 from Fig. 3b, whereas for a capacitative one,6 = 1. 


We can now calculate the steel volume and the number of turns on the a-c winding: 


a | 
= ppt — 108, cm’, 
ee (37) 
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st 
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U, 108 


" 4fSs, Brin 


Ww. (38) 


where Bin ts in v-sec /m* , H.wmax in amp-turns/m, and S (core section) in cm*, 




















Fig. 15 


In conclusion,we observe that as the operation is highly reliable the dimensions are somewhat overesti - 
mated, so the method cannot be recommended for high-power amplifiers, particularly if transformer steel is 
used. 


The resistance of the rectifier can be allowed for approximately from its voltage-current characteristics, 
The resistance calculated for the maximum load current should be added to Rj. A more precise method is 
given in[1]. The resistance of the windings can be added to R, in the same way. 
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CONTACTLESS DETERMINATION OF SPECIFIC ELECTRICAL RESISTANCE 


G. G. larmol'‘chuk 


The utilization of a high-frequency current when determin- 
ing the specific electrical resistance of matter is considered. The 
circuit of a bridge for the contactless determination of specific 
electrical resistance is described,and formulae to determine the 
sensitivity of the bridge circuit in two cases of measuring are given. 


1. Physical Principles of Methods 





The contactless method of measuring the specific electrical resistance p , developed at the Institute of 
Automation and Remote Control of the Academy of Sciences of the USSR, is based on the utilization of eddy 
currents which arise in a conducting material when the latter is placed in an alternating magnetic field. 


The alternating magnetic field is set up by an induction coil which is supplied with a high frequency 
current. A body placed in the magnetic field of the coil changes the latter's parameters. This change in the 
coil’s parameters depends on the geometrical dimensions of the body and its specific resistance p. The coil's 
parameters may be measured by any known method. 


Let us examine this phenomenon more closely, Figure 1 shows a uniform cylinder of a meterial having 
a specific resistance p , lying in a high frequency magnetic field, Circulating currents are induced in the body 
of the cylinder in planes perpendicular to the vector of the magnetic field H. 


Due to the skin effect, the current distributes itself nonuniformly through the cross section of the speci- 
men. The greatest current density is found on the surface of the cylinder. The distribution of the current 
density through the cross section of the test piece is shown in Fig. 2, where dg is the surface current density, 
8 is the skin depth, $ D is the radius of the cylinder. 


Since the use of such a curve is difficult, the 
effective action of the current in the body of the cylin- 
der is replaced by the current density q9, spreading in 
the cylinder to a depth 6. 


In accordance with known data [1-3], the skin 


depth is 


where p is the specific electrical resistance of the 

Fig. 1 Fig. 2 material in ohm mm*/m, y is the magnetic permea- 
bility of the material, f is the current frequency in 

cps. 
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For nonmagnetic materials (of which more will be said later) y = 1 and formula (1) takes on the form: 


o= 903 V3 MM, (2) 


The current induced in the body of the cylinder sets up its own magnetic flux, which is determined by 
the skin depth6 and the dimensions of the cylinder. Hence the cylinder may be considered to be hollow, 
having a specific resistance p and wall thickness 5. Such a body will have its own resistance R, and inductive 
reactance wh, where Lg is the inductance of the cylinder. In this case the cylinder may be represented by an 
equivalent circuit, consisting of series connected Rg and Ls (Fig, 3). 


Such a representation of the phenomenon 
occuring in the body of the cylinder reduces it 
to a case which is studied in the theory of 
coupled electrical circuits, In this case, the 
cylinder and the induction coil, in whose field 


po the cylinder is located, are effectively, in- 
‘ : é t, 3m Bu, # Ry 
2 6 * . ductively coupled circuits. 
tt 


Fig. 3 Fig. 4 Figure 4 shows the circuit diagram of 

inductively coupled circuits. Circuit 1 con- 
sists of a coil, having a inductance Ly, con- 
nected in series with a resistance R, and capacitance C,, where E is the current source. Circuit 2 represents 
the equivalent circuit of the cylinder, where M is the mutual inductance between the two circuits. 








From theory [4] is known that, in the presence of the mutual inductance M, the effect of the second cir- 
cuit on the first is to change the resistive and reactive components of the first circuit. 


The effective impedance of the equivalent circuit, i.e., the impedance with the effect of circuit 2 on 
circuit 1 taken into account, will be 





2 M2R ‘ wo'*M2y 
Zio = Ry + Fo + 7 (X1 —S) 
1o0e i+ | Z|? +] 1 |Z, |? , (3) 
where Z, is the effective impedance of circuit 2, Z, = Rg + jXg, Xy = why, Xg = uly. 
In Equation (3) the second term on the right side represents the magnitude of the added resistance. 
A R = o'*M2R, aie w?M?R, 
Tea RE (Lay (4) 


If the skin depth 6 is small compared with the dimensions of the cylinder, then the inductance of such 
a cylinder will be affected very slightly by the magnitude of p, which affects the skin depth, while the re- ni- 
sistance Rg will be proportional to vp. 





In fact, Rg, for 5 << , may be calculated from the formula 


D 
Ry = p--- (5) 


Substituting the expression for 6 from Equation (2) into (5), we obtain 





ree exD V 7 PV of 
zs. 503hV p 107 5.03 x 10% * (6) 
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Here D and h are expressed in millimeters, P = 1D is the circumference of the cylinder. 


Thus the magnitude of the resistance R, of the cylinder will depend on the specific resistance of the 
cylinder's material, which in its turn will bring about a change in the magnitude AR,, i.e., the transfered 
resistance in circuit 1. 


To find the conditions for a maximum AR,, let us differentiate (4) with respect to R, and equate it to 
zero: 


d(AR) _ _OMP(oLsY—R __ 
ah; [A2 + (wL,)*]? (1) 





Thus we find that AR, is a maximum when Rg, = ul». 
It follows that AR, is a function of the expression Rg /ul,. 
Let us set Ry /wh, = n; then Ry = nuly. 
Substituting this expression for R, in formula (4): 


_ @M?wLon 
AR, = (oLa)® (a? 4-1) * 





and thus 


AR,L. n 


oM? ~ nt +T° (8) 





Figure 5 shows the dependence of rie on n. It is apparent from the curve that the maximum occurs at 


n= Lé 


Equation (4), when AR, is a maximum for a given 
w, M, and Lg, may be rewritten as: 


M2 
Rema o 





The more the R, will differ in the sample from 
the inductive reactance uwl2, the more the transferred 
resistance AR, will differ from AR,max. 





G0 Vi / 1 0n 
In the measurement of p for samples having equal 
dimensions, and considering wl, constant at least with- 
in range of the measured p's, a single valued solution 
requires that either the left or the right branch of the curve in Fig. 5 be utilized. From the point of view of 
greater sensitivity it is recommended that either region A or region B be used (Fig. 5). 


The point of inflection occurs where R, = wh». Substituting the expression for Rg from Equation (6), we 
obtain 


PV of 
“6,03 x 10°h = wl. (10) 
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From Equation (10) it is possible to obtain the value of Por, 1.e., the value of p corresponding to the point 
of inflection in the curve of the transferred resistances: 


h \2 
Per = L(+) (10 (11) 


From Equation (11) it is possible to determine that critical frequency, for a given value of p , which cor- 
responds to the point of inflection in the curve: 


Jer “ie (x) 10~™ cps. 


4 


(12) 


In Equations (11) and (12) Lg is expressed in henries , P and h in mm, f in cps, and p in ohm mm*/m. 


The nature of the dependence of AR, on h was verified experimentally. Measurements were carried out 
on a QO-meter. 


The experimental curve for a prismatic sample, measuring 5 x 6 x 24 mm, at a frequency of 150 kc, is 
given in Fig. 6. 


The dependance of the position of the point of inflection,in the curve of the transferred resistance,on fre- 
quency can be seen in Fig. 7, where values of the figure of merit Q = 1/ wCR are plotted along the ordinate. 
As R increases, Q decreases. Hence the maximum transferred resistance corresponds, at a given frequency, to 
the lowest value of Q. The curves are given for prismatic specimens measuring 5 x 6 x 24 mm, for different 
values of p and for different frequencies. It is apparent from the curves that, with increase in frequency, the 
minimum Q moves in the direction of larger p , which agrees with the condition for obtaining ARymax (Rz = wl). 


From the curves it may be also seen that for a sample measuring 16 x 20 x 50 mm, at a frequency of 
1 Mc the minimum Q occurs at a much larger value of p than for a sample measuring 5 x 6 x 24mm. 


From what has been said, we can establish rules 
which must be followed to obtain single-valued mea- 


3 





If the rising portion of the curve in Fig. 5 (region 
A) is used in the measurements, then the maximum value 
2 ? aia we Pmax, for a given frequency, must be less than the value 
of Por, For the dropping part of the curve (region B), we 
must accordingly have Pmin > Pcr- 


























’ Ly X\ Insofar as, with increase in frequency, the point of 
pete inflection in the curve [cf.Equation (12)}moves in the di- 
rection of larger p , then the current frequencies must be so 
40s at 7 0 on chosen that the point of inflection of the curve falls with- 
Fig. 6 in a controllable region, enclosed by Prin and Pmax- 


Using formula (12) we find the frequencies fpryjn and 
fPmax, at which the point of inflection of the curve will 
‘correspond tO Prin aNd Pmax- 


For measuring purposes, any frequencies may be used (from a single-valuedness point of view) except those 


lying between fpmin and fpmax. Hence for region A the frequency must be greater than fpmax, while for region 
B—less than fpmin. 


Thus the conditions for single-~valuedness require that the following restrictions be observed in measuring 
p: for region A (Fig. 5) Pmax < Pcr and f >f Ppyax, for region B— Pmin> Per and f<Pmin- 
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2. Methods of Measuring the Specific Resistance 





The physical basis of the contactless method of measuring p which was discussed above, shows that the 
problem is reduced to the measurement of the transferred resistance in the coil circuit. This resistance may 
be measured by any of the well-known methods; in particular, bridge circuits may be used. 


The bridge circuit shown in Fig. 8 was used at the Institute of Automation and Remote Control to check 
the specific resistance of small samples [5] using the contactless method. The following symbols are introduced 
in this diagram; 


Lj is the variable inductance in balancing arm of 
the bridge, Lg is the measuring coil having inductance 
Le, ty and rg are variable resistances, C, and Cy are vari- 
able air-gap capacitors, Tp is a coreless transformer for 
isolating the bridge circuit from the VTVM circuit, KV 
is a vacuum tube voltmeter, Up is an applied source 
voltage of 7-8 volts at a frequency f = 1 Mc, Ly and L, 
are filament lamps. The balance conditions for the bridge 
are; 


a (13) 


From the balance conditions we can write 


LC, = LC. (14) 


Equality (14) is independent of the frequency w. 


If we pick a value of w such that wLjC, = WLC, = 1, 
then conditions of series resonance will exist in both arms 
of the bridge. 


The use of the resonance effect is convenient to the 
extent that the increased sensitivity which accompanies it 
permits the use of a VTVM as an indicator. 


As practice has shown, an applied potential Up = Tv 
and f = 1Mc can produce out of balance voltages of tens 
of volts. 





Two methods of measurement may be used with the 
given bridge circuit. 


Fig. 7 


1. By changing the parameters of the bridge arms, 
resonance is obtained in both arms of the bridge (this is 
determined from the filament lamps L; and Lg) at the same time as the bridge is balanced (indicated by the VTVM). 
Then a master sample, having a known resistance , corresponding to one of the limits of the p being measured, is 
introduced into measuring coil L3. The deflection on the VTVM is noted. Next,the sample in coil Lj is replaced 
by another master sarnple having the same dimensions but having a specific resistance corresponding to the second 
limit of the p being measured. The new voltmeter deflection is noted. Following this, all samples having the 
same dimensions as the two master samples have a p obtained from the formula; 


_ U,~-0, 
Px = 7, (Pa — Pi) + Pr» (15) 


where U, and p, are respectively, the out of balance voltage corresponding to first master sample, and the latter's 
specific resistance, U, and pz, — the same for the second master sample, Ux and py — the same for the test sample. 
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Formula (15) is derived on the assumption that U is a linear function of p in the established limits of 
measurements, and that this holds true for samples which, when placed in the test circuit, produce voltmeter 
deflections greater than those corresponding to the master samples. 


2. For the initial calibration of the bridge a master 
sample is chosen which has a specific resistance p, lying 
beyond the limits of the p expected for the samples to be 
measured, In this case either py > Prjgys OF Py < Prin 
where Prin and Prax are the assigned limits of mea- 
surement. 





Next,sample py is replaced by another master sample 
with specific resistance py. The value py is chosen so that 
the range of measured pis bounded by p, and pz, i.e., if 
P2 < Pmin. then Pg > Prax, but if py > Prax, then 
P2 < Prin (Fig. 9). 











When the master sample having pz, is introduced, 
the bridge is unbalanced and a voltmeter reading is 
registered. Following this, the resistance of every un- 
known sample p, can be determined from the formula 





U z 
me x ee ohm -mm 
Px = Pi + Us (P2 Px) m ’ (16) 


if pp >p,, and 


x 


17 
Px = Pr — 7 (Pi — Pra) a 





if pp < py. 


Formulas (16) and (17) are derived on the assumption that the voltmeter readings are proportional to the 
difference p, — pz. 
The sensitivity of the bridge is greater in the second method than in the first. 


As was pointed out above, filament lamps L; and L, are included in the arms containing the inductances 
Lj and L3 to facilitate tuning the arms to resonance. At resonance the lamps glow. The sensitivity may be in- 
creased by shorting the lamps with keys (after resonance is achieved) and thus decrease the resistance in arms 
1 and 2. 


3. Sensitivity of the Bridge 





When both arms of the bridge are tuned to resonance, and if conditions (13) are met, then the bridge will 
be balanced. Let us pick C, = Cy, 


If the second method of measurement is used, i.e., if the bridge is calibrated using the master sarnple p,, 
the bridge will be put out of balance when the master sample is replaced by a test sample having a different 
p. In the absence of current in the indicator loop (diagonal CD in Fig. 8) the voltage across capacitor C, will 
be equal to (at voltage resonance); 





ry joCs ’ (18) 
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where r, is the total resistance in the arm ADB, Up is the voltage delivered by source of potential, Cy is the 
value of capacitance. Analogously for branch ACB 





U 
en et Paint (19) 


re jo, ¥ 
where ry is the resistance of arm ACB. 


The voltage Ug across the detector diagonal of the bridge will be equal to the difference between the 
voltage drops (ignoring the effect of the current in the detector diagonal): 


4 _ ~~ P \rijaCs rea, }* 


Since Cy = C, by choice, then 


U 
wpe (tt). 
Jas \ ry rs (20) 
When the master sample in coil Lj is replaced by a sample with a different value of p , the resistance 
in arm ABC will become equal to rg = rg + Arg, where Ar, is the change in the transferred resistance in the 
circuit of Lg due to change of sample. 


In this case expression (20) takes on the form: 


Ve~ Fate (EER) 


When the bridge is balanced with the master sample in the circuit of Lg, the equality rg = r, exists and 
one can write 





em oh, 4 |. 


(71 + Ars) (21) 
he brid nbs 
The ge sensitivity is S = Oar, 
Taking the derivative of (21) with respect to Ars, we find 
= U 
Rept siglmallees. tel 
. jaCs (r,-+ Ar,)* ° (22) 


If ry >> Ar, then we can assume that S does not depend on Ar, and remains constant, while the voltage 
Ug, according to equation (21), will be proportional to Arp. The scale will be linear, Usually the sensitivity 
increases with decrease in 1, i.e., with increase of Q-factor of the induction coils Lf and Lg. In fact, taking 
the modulus of formula (22) we get 











ral U 1 UpQ 
Csr; A 2 Ar, \2’ 
n(i4tt) on (1+) ie 


where Q, = 1/wCgr, is the Q factor of circuit ADB. 
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To determine the dependence of S on p we sub- 
stitute (6) in (4); 








Poe tae 
ANT NJ = 5.03 x 10% 
| » Ar, => P V of 2 7 - 
X Y sn x a + (ols? ae 
q 
MNS 
Piz Prin = =Pmax Fas The calculations may be simplified by determin- 
Fig. 9 ing the bridge sensitivity S under the condition of maxi- 


mum transferred resistance, i.e., under condition (7). In 


such a case Equation (27) can be rewritten in the follow - 
ing form: 


2172 
Ar, = «*?M25.03 x _10*h 
2P V of 





(25) 


When a sample having a specific resistance p is replaced in the measuring coil by another sample with 


specific resistance p + Ap, the change in the circuit's resistance becomes equal to the difference between the 
transferred resistances; 




















Ar. ex 2:08. 10°o*M*h (fase 1 ) 
: 2PV 7 Ve Voe-+Ap: (26) 
Th i : : in (26) i imately equal t 
e expression n $ approximate e O: 
4 i; Fee ee ' 
4 1 i il 
Veo Vetde_ »Ve (27) 
Let us introduce the symbol 
5.03 « 105a*M*h 
A= — 
4P Vi (28) 
Then we will have 
, Ape 
Ar, = A———. 
; eVe (29) 


We substitute the value of Arj from (29) into (21) to express the potential across the measuring diagonal 
of the bridge: 


Ty ae Ae ; 
g 700s ry (rp Vp + Adp) (30) 





257 








To find the modulus of the bridge sensitivity as a function of 4p, we take the derivative of expression 
(30): 


Up5.03 xX 10a*M*%hrip Vp 4PV 7 


Sp = = = = . , 
; wary (4P V {rip Vp +5.03 X10°w*M*hAp)# (31) 


0A 





la" 





If, in the denominator, the second term in brackets is smaller than the first term by an order of one or 
two (as a result of the smallness of M”), then 


: pee Upn V7 M*h5.03x 105 
2Cr? Po Vp a 5 (31°) 


Equation (31°) can be rewritten in a different form by introducing the figure of merit Q, = 1/ uCgry: 








5 Uphr V 7@Q,5.03 x 105M? (32) 
R®: 2r, Poe Vo ; 
PVvpi 
Maximum sensitivity occurs when R, = wl». Replacing in Equation (32) 503x 10h R, by ul, and 
M2 by "LiL, we obtain ’ 
Ss Up Kel, Q: (33) 


°max v's 4rie 


It follows from Equation (33) that to increase the sensitivity it is profitable to increase the inductance of 
the measuring coil, the Q-factor and the coupling coefficient. The sensitivityincreases as well with decrease 
in Ty. 


All the conclusions are arrived at on the assumption that the distributed capacitance of the inductive coils 
is small and can be ignored when compared to the capacitance Cs. 


In practice, the bridge sensitivity S, varies by fractions of a volt to several volts per unit of specific re- 
sistance, depending on actual conditions. Such a sensitivity permits the use of a VTVM as an indicator without 
an intermediate amplifier. 


The transformer (cf.Fig. 8) is used in the circuit mainly as an isolating element. However, the bridge sen- 
sitivity may be increased by increasing the transformation ratio of the output transformer. 


4. Range of Application of the Described Method [9-12] 





The described method may be used to measure small samples of nonmagnetic materials having a volume 
of 20-30 mm’, The measurement of p for such samples by usual methods using a dc source, presents great 
difficulties, 


The method may be used as the basis for the development of apparatus for the determination of the p of 
materials in ingots, manufactured objects, etc. It is also useful in the design of automatic sorters of items by 
their specific resistance. Using phase discriminators one can control the dimensions of the items, determine the 
presence of logitudinal cracks appearing on the surface of the items, and thereby control the quality of production, 
etc, 


It is possible to devise an absolute method of determining p by measuring the transferred AR, and AL, 
once a method is developed for the accurate determination of the magnitudes R, and Ly in relation to the fre- 
quency and the dimensions of the samples, 
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5. An Apparatus Using the Contactless Method for the Sorting of Carbon Brushes 
by Their Specific Resistance 








The specific resistance of electric carbon brushes according to GOST 2332-43 varies between 0.05 and 
60 ohm mm*/m. 


The resistance of the brushes to current Rgpch depends on p and their dimensions. It also varies over a 
wide range. Simple calculations yield values of Rshch for brushes in the range from 1x 10° to1 x 10°" ohms. 


Due to the small dimensions and resistances of the brushes, the measurement of the specific resistance of 
carbon brushes represents a difficult problem. 


The specific resistance of brushes is measured either by the method of separated circuits or by a bridge 
method, However, the variation of the current density along the cross section of the brush, as a result of the 
effect of the quality of contact between the brush and the connecting plates, which arises in the method of 
separated circuits, or else the contact resistance between the current.carrying plates and the brush, introduce 
significant variations which affect the magnitude of the experiment errors in the method, Practice shows that 
these errors may sometimes reach values of 100% of the resistance being measured, 


Neither is the problem of control in mass production solved by the rho-apparatus developed at a branch 
of the research institute for electric cells and electrocarbon, which is based on the method of divided circuits 


[7]. 


These difficulties in measurement (control) were overcome by using the above-described contactless 
method of measuring the specific resistance. Models of the bridge described in part 2 were built for this pur- 
pose. A photograph of one of these units is shown in Fig. 10. The measuring coil into which the brush being 
measured is placed can be seen on the right. A single unit comprises the bridge and a single-tube high fre- 
quency signal generator. Power is derived from a regulated rectifier which is fed from 220 volt mains. The 
signal generator frequency is chosen between 0.6 and 1 Mc for different models of the apparatus. 


Results of measurements made with the different models were found to be easily reproducible and the 
sorting speeds were high. 

















Fig. 10 


The inadequacies of the bridge include the necessity for a preliminary warm-up before measurement, and 
the necessity of determining the specific resistance of the master brushes in a laboratory by an absolute method 
(amalgam or rho-apparatus) for each type of brush being controlled. 


The results of control-measurement on the given equipment sometimes diverge from the results obtained 
using de current methods owing to different laws governing the evolution of the mean measured value, However, 
the results are much more reproducible on this equipment than on old methods, 
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The apparatus is used to control stocks of brushes prior to their copper plating. However, on principle, 
the method may be applied to the control of copper plated brushes, without ruining them. 


The equipment will work either if the brush is introduced into the measuring coil or if it is placed on 
the face of the coil, thus facilitating the design of an automatic sorter. 


When brushes are produced to the third grade of accuracy, the theoretical rejection error is approximately 
8%, of which 5% is due to the effect of the tolerance in the length of the brush. 


The circuit of the apparatus is simple, and it is suitable for manufacture in a factory KIP. The most dif- 
ficult part of the work involved in assembling the bridge is the initial tuning of the high-frequency bridge. 
Equipment of this type has been used successfully for over a year and a half at one factory producing carbon 
brushes. 

SUMMARY 


The theoretical conditions which were discussed in the paper point up the possibility of utilizing, for con- 
trol purposes, currents, which make it unnecessary to have contacts between the object being measured and 
the elements of the measuring circuit. The developed models of a device for the sorting of objects by their 
specific resistance have shown their advantages over sorting equipment using a dc current. 


This method may be used in production for the continuous control of some parameters which determine 
the quality of the product; the specific resistance, diameters of objects, the presence of longitudinal cracks, the 
wall thickness of thin-walled cylinders, etc. 


This method can be recommended for the development of automatic control and sorting equipment, and 
in some cases it may serve as the basis for developing pickups for application in automatic regulating systems. 
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A LINEAR INDUCTION POTENTIOMETER FOR INDUSTRIAL USE 


O.I1. Aven, S. M. Domanitskii and Io. M, Pul'‘er 
(Moscow) 


A description is given of an induction potentiometer with 
a linear dependence between the output voltage and the shaft 
rotation in a range close to + 90°. A method of computation is 
given, and also an experimental method for investigating similar 
potentiometers. Experimental results are provided. The potentio- 
meter described could serve as a standard for industrial types. 


INTRODUCTION 


The slide-wires which find such broad use in industrial automatic devices have essential inadequacies, 
The presence of sliding contacts, which contribute to the unreliability of the slide-wire functioning, also affect 
the reliability of the apparatus of which the slide-wire is a component. In spite of its apparent simplicity, the 
accurate preparation of slide-wires is quite complex. The use of slide-wires together with contact-less apparatus 
in, for example, the positioning feedback loop of a contact-less executive mechanism [1], leads in fact to a 
de-emphasis of the contact-less parts of the devices. The inadequacies of slide-wires necessitated the search for 
a means of constructing contact-less devices which still possessed characteristics analogous to those of slide-wires, 


In contemporary apparatus destined for use in industrial automatic control systems, there are sometimes 
used, instead of slide-wires, selsyns or inductive transducers. However, selsyns have a maximum rotational 
angle, corresponding to the linear portion of their characteristic, of about + 15-20°, which can only be extended 
to + 50-60° by adding special windings, which frequently leads to all sorts of troubles, Small-size inductive 
transducers are usually manufactured with forward motion and are designed for small disturbances. Moreover, 
there is a transformation of the forward motion to a rotary one which complicates the kinematics of the appara- 
tus. Thus, there are no contemporary contact-less variable transducers of the rotary type for industrial use which 
possess large rotational angles (180-240°). 


The present paper describes an induction potentiometer designed, in particular, to replace slide-wires in 
control and regulation systems which have a maximum rotational angle of 180° (490°). It should be mentioned 
that the inductive system on which this potentiometer is based allows a maximum angle of 240° (4 120°) to be 
obtained, The problem in developing this potentiometer was to make it universal, suitable for use in secondary 
tools, regulators and executive mechnisms. 


With this problem for a starting point, the following special requirements were formulated for the potentio- 
meter. 


1. The potentiometer should function in equalizing designs, similar to the schematic of Fig. 1. 
2. The potentiometer must be designed for an ac source voltage of 50 hertz frequency. 


3. The potentiometer 's output voltage must correspond to the input parameters of the standard industrial 
control system amplifiers. 


4. For rotational angles close to 180° (+ 90°) or 240° (+ 120°), the error in amplitude must not exceed 0.5%. 





5. The phase angle between the input and output voltages must not depend on the angle of potentiometer 
shaft rotation. 


6, In construction and dimensions, the potentiometers should be compatible with control and regulation 
devices (executive mechanisms, regulators, secondary tools, etc.). 


Potentiometer Construction 





The induction system described in [2] was taken 
as the basis for the functioning of our industrial linear 
induction potentiometer. 


Figure 2 gives the physical design of the linear 
induction potentiometer with a + 90° range of shaft 











yee, | rotation. On the same figure is provided the design of 
Ly, mi the magnetic link, with the principal dimensions noted, 
Fig. 1 and an indication of the fluxes, ®, and , for different 


magnetic flux paths, as determined by the given position- 
ing of the windings and the choice of the current direc- 
tion. 


There is a pair of windings W., and Wg, on the potentiometer rotor, and an analogous pair Wy, and Wy, 
on the stator. 


The primary windings W,, and Wy, comprising the exciting winding, upon rotation of the rotor in the 
range + 90°, change their orientation with respect to windings Wg, and Wg, comprising the secondary winding. 
At this, the flux in the secondary winding changes from zero to a maximum value, at some phase angle. The 
secondary winding is loaded by load resistor z,, . 
































Fig. 2. Physical design of the potentiometer. 


The construction of the potentiometer is clarified by the overall view (Fig. 3). In the slots of the non- 
magnetic housing of the stator 1. are embedded the blocks at the back of the stator 2 composed of drawn 
sheets of electrical steel. These sheets are bonded to CF. The layer of adhesive serves simultaneously as in- 
sulation between the laminas. Also inserted in the stator housing are the lateral 3 and intermediate 4, disc 
block which, together with the back, make up the Ill-shaped magnetic stator. The stator winding 5 is wound 
on a former and inserted in the stator assembly, 









































Fig. 3. Overall view of the potentiometer. 


The rotor is assembled the same way as the stator and, analogously to the stator, has the III-shaped mag- 
netic characteristic (8, 6 and 7 are, respectively, the rotor back, side and middle disc blocks). The rotor wind- 
ing 9 is also wound on a former and inserted in the rotor assembly, the winding endings being led along the 
uncovered shaft, The potentiometer construction is intended to allow the mechanization of the winding on the 
stator and rotor. 


It should be mentioned that, as shown in [2], there is the possibility of so laying out the primary and sec- 
ondary windings on the stator that one may obtain a completely contact-less system for arbitrary angle of ro- 
tation. However, in a system with windings of the rotor, with rotation angles less than 360°, it is possible to 
have the winding endings made of flexible conductor, thus securing a contact-less systern. The construction of 
a potentiometer with windings on both rotor and stator, in comparison with a design having windings only on 
the stator, possesses a host of advantages — effective utilization of copper and steel by virtue of the complete 
encircling of the primary and secondary windings along their entire length. This latter fact accounts for a 
minimum of end effect, thus increasing the effective angle of rotation, and permitting a high degree of magnetic 
symmetry in the system. This in turn effects a much higher accuracy of potentiometer behavior. 


Method of Potentiometer Calculation 





The method for carrying out the engineering calculations for the potentiometer is based on the given 
technical requirements. We make use of the basic relationships given in [2]. Below is given the method for 
computing an induction potentiometer for the case where one is given the supply voltage U, the maximum 
amplitude of the secondary voltage Ugma x the load resistance zg and (+ 5 Ug)max the active component 
of the maximum potentiometer error (coinciding in phase with U,), taken as the ratio of Ugmax with respect 
to unity. 


From the results of the computations are determined the physical dimensions of the potentiometer, the 
number of turns and the conductor diameter for the primary and secondary windings, The magnitude of the 
primary and secondary currents is found, as well as the output voltage when the potentiometer is loaded. 


The computations are carried out in the following sequence. 


1. We find the ratio of the number of turns, k = wg /w,, of the secondary and primary potentiometer 
windings and also the inductive reluctance xX». Using Relationships (1) and (2), which appear in [2], 
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(1) 


(+ 8Us)max = 0.25Re ker 


rT 
“omax_ __k 


/_ ma" (2) 





we obtain the system of two complex equations (3) and (4), with repsect to the argument of Ugmax (the angle 
3) and the coefficient k; 


Top + 104 





U2 max(COS 2 + J sin o,) = U,k 


(8U2)max (734 + 23) 


I 


y = k*z, = 2.6 





(4) 


For deriving Equations (3) and (4) from Relationships (1) and (2), we take r, = rg = 0, Xy5 © 0.3Xq and 
Xgs © 0.3x¢k?, Moreover in the expression for A(1) in Equation (1), we take zy = 0. Then 











ai)= (22 +1) (ge41)4 Mats, 


£9) 


The following notation was employed in giving the previous equations; 2,5 = 2 (r, + jxys) is the total re- 
actance of the primary winding (stator), Z, = 2 (tg + jxgs) is the total reactance of the secondary winding (rotor), 
Zg, = Tel +jxXg_ is the total load reactance, r, and xs are the active and inductive reactance of one coil of the 
stator winding, rp and Xg; are the active and inductive reactance of one coil of the rotor winding, rg, and xq, 
are the active and reactive load reactance, and w, and wy are the number of turns in the stator and rotor wind- 
ings. 


From Equations (3) and (4) we obtain 











ke = Z2max Vs Bry + (1-32, 1 + 1. 38yP 


sy rit ty - 
Here, y is immediately determined for (4) and, in its turn, 
a yt as 
Zo — k2 . (6) 


It should be mentioned that if zy = oo, or ry, then Equation (1) reduces identically to zero. In this case, 
one has only Equation (2) for a starting point, together with (3) and (5), which are derived from it, from which 
it is clear that, for y = 0, the value of the reactive reluctance can be arbitrary. As will be shown below, X9 
in this case is determined from practical rational relationships deriving from the geometry of the induction 
system and the given current requirements. 


2. We find the current of the unloaded potentiometer; 


(7) 
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where OL AL 2169" 





4ui= 2171 + 7 (2s + z,)], 711 = 0, Ly, = 0,32, 


3. We assign the maximum induction, Bmax, which is obtained in the magnetizable disc in the section 
immediately adjacent to the back of the rotor. For transformer steel of type E42 (E4AA), which will essentially 
render neglible the influence of nonlinearity of the steel, we take Byyax = 2000— 3000 gauss. 


4, We determine the magnetizing force necessary to generate the induction Byyax, not taking into account 


the reactance of the steel sections: 


F, = Ly ?1 = 0.8 x 2Bmax(“) é. 


d (8) 


The magnitude of the air gap, 5 = 0,2-0.3 mm, and the ratio, d,; /d = 0.6-0.8, were chosen from con- 
siderations of the physical configuration (the dimensions d, and d are the diameters of the stator bore, as shown 
in Fig. 2). From Equation (8) we determine the number of turns in the stator winding: wy = F,/Iyui- 


5. We find the magnetic flux connected with the stator winding: 


1 4.44fu, 


(9) 
Here, Us = IULXo = LyL“AuL/2-6 = U;/2.6. 


6. We determine the width of the blocks ag and a, (Fig. 2), of the stator and rotor magnetic conductors, 
as a function of the diameter d, given by physical considerations; 


Bryaxth, \di/’ ~*~ 2 (10) 


(it is desirable that a, > 105). 


Here, r = 1d/2,ke » 0.9 is the charging coefficient of the steel blocks. 


1. We compute the cross-sectional area of the backs of the stator and rotor magnetic conductors; 


gy. eee. 
Bmax (11) 
from which the basis dimensions of the backs may be found. 


8, We determine the area S, = Iy,yy /A@nd the diameter, d,) of the stator winding conductors by start- 
ing from the given current density, 4. 


9. We find the beight h of the aperture (Fig. 2) by taking the width of the aperture to be c ~ 2a; 


s 
Pa (12) 
c 2k a1 





Here, Sap is the aperture area, k, ~ 0.3-0.35 is the charging coefficient of the aperture copper. 


10. Knowing all the magnetic conductor dimensions we determine, from known formulae, the active 
reactance of one stator coil. 


11. We find the number of turns of the rotor winding, wz = kw,, the area of conductor S, and the react- 
ance of one rotor coil fp. 


12. After the determination of the winding data, the current and voltage in various regimens, taking 


into account the ohmic resistance of the windings and the influence of the load, are found from Formulae (13), 
(14) and (15) given in [2]. 
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The relations referred to will be complex, thus allowing the determination both of the real values of 
the current and the voltage, and also their phase displacement. 


The current in the stator winding 
h« iy 


; 2) 2ok* a i’ 
244 + 2% 1 wits deck oo | (13) 








where a / Ompgx is with respect to the angle of rotation, 
Zq = 2y5 + 2g) = 2 (72 + fps) + 294 
The current in the rotor winding 


i 2jxokl, a 
2 ‘ N 
22 + 2/2ok, a 





: (14) 
max 


The magnitude of the stator and rotor winding currents without loading, and short-circuited, may be found 
from Expressions (13) and (14) if one takes, respectively, Zg= co and Zg = Zgs, 


The voltage across the load ae 


. a ° 
U, = — = Usk. 


205 z z,,*? 2) xyk* a 
pa ae 18 pau 
bs * (zs +1)+ 2951 t 2] |! (axa) | sine 














Here, a is the angular rotation of the potentiometer rotor, a4, is the maximum angular rotation of the 
potentiometer rotor. 


The magnitude of k, and the active part of the maximum error ( 5Us)max may also be specified from 
Expression (5) and (1). 


In the cases where Zg, = © Or Zg, = rg, the computation is carried out in the following manner. The 
magnitude of the required current I,,jj, and the diameter d are given from practical considerations once Uj, 
Ugmax 4nd Zg, are known. The number of turns w,; may be found from Equation (8), and the reluctance X, 
from the expression, x9 = U, /2.6l,j,. Further computations are carried out in the same order as above, taking 
into account the rational relationships in the geometry of the potentiometer, a; >105, c/a; #2, ag /a, ~2, 

d, /d = 0.6-0.8. 


EXPERIMENTAL METHODS 


Below are given the results of experimental investigation of the linear induction potentiometer. The 
potentiometers were calculated on the following data; U, = 80 volts, Ugrjax = 7 volts. They had the following 
dimensions; D = 70 mm, d = 50 mm, a; = 3 mm, ag= 4mm, 6 = 0.3 mm. Photographs of a disassembled 
potentiometer are given in Fig. 4. 


The linear induction potentiometers were investigated both as individual components, and as functioning 
parts of equalizing designs. As shown by experiment, the empirical curves, Ip = f (a), sc = f (a) and Igsc = 
f (a), for individual potentiometers, have the same character as the computed curves constructed from Equations 
(13) and (14). The maximum discrepancy for the currents, I,1jj, and Iysc,, nowhere exceeded 10-12% and for 
the current Igg-,, nowhere exceeded 5-7%. Satisfactory agreement with computed values was also given by the 
empirical curves for the phase displacements, gj, = f (%), gay, = £(&), Yisc = f(a) and gegc = f (a), 
where yyy, is the angular displacement between the primary current Lu. and voltage, U,, vectors, ygyy is 
the angle between vectors 0, and Usui, vic 1s the angle between U, and isc and ¢gsc is the angle between 
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U, and iggc. The angles yy, and yy do not, in fact, depend on the position of the potentiometer rotor. 
The magnitude of the angles was measured by means of a Ts50 vectormeter. The maximum deviation between 
the experimental and the computed values was 5° for gp and gyygc and 12-15" for gg and ¢ggc. 























Fig. 4. 





Fig. 5 


Figure 5 shows the dependence of the secondary voltage on the angle of rotor rotation for one potentio- 
meter, with load resistance Ry, = 500 ohms (curve 1) and R; = (curve 2). The minus sign for voltage Us is 
conventionally chosen so as to show phase variation of the output voltage from 180°. 


Figure 6 gives the initial portion of the curve for R, = @. The angle of rotation was established with a 
high degree of accuracy by means of an optical graduating device. The minimum “residual” voltage, Ugmin, 
corresponding to the angle a = 0, for different potentiometers was measured at from 12 to 35 millivolts, 
comprising 0.2-0.5% of Ugmax. It is necessary to mention that the residual voltage contained fundamentally 
the higher harmonics, The relationship between the voltage at the fundamental harmonic and the voltage at 
the higher harmonics in the curve of the secondary voltage, Ugmin, depends only on the quality of the mag- 
netic material, and on the care with which the potentiometers are prepared, The composition and magnitude 
of the higher harmonics depends, in significant measure, particularly on the quality of heat treatment of the 
steel plates of the magnetic conductors, but the magnitude of the residual voltage at the fundamental harmonic 
also depends on the accuracy of the preparation and assembly of the potentiometers. Figure 7 gives the depen- 
dence of Ugmin, Ugmax 29d Ugmin /Usmax for one of the potentiometers on the supply voltage U,. As is 
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clear from Fig. 7, for small and large values of supply voltage, the quantity Ugmin/Ugmax increases, This 
stems from the fact that, for weak and for strong fields, the influence of the nonlinearity in the steel character- 
istic is greater, thanks to which the difference in e.m.f.'s of the potentiometer rotor winding coils is increased. 
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For various regulators, particular interest inheres in the magnitude of the “residual” unequalized voltage, 
corresponding to the actual adjusted state of the potentiometer rotor in an equalizing circuit, This voltage is 


fed to the input of the regulator's electronic amplifier, and should be as low as possible in order to minimize 
the static error of the regulating system. 


V-esidual » millivolts 
100 > 














Fig. 8 


Figures 8 and 9 give the curves Ures = f (a), taken from equalizing circuits for U, = 80 volts and Ry = o. 
In the first case, the stator windings of two potentiometers were connected in parallel, The rotor windings were 
loaded with a resistance, R, = 10 kilohms. In order that the rotor winding voltages coincide in phase a resistor 
R, = 140 ohms was connected in series with the stator winding of one of the potentiometers. For the series-con- 
nected stator windings (Fig. 9), one of its windings, with the same purpose, was shunted by a resistor Rg = 330 
kilohms, 


In the potentiometers under consideration, the output voltage U, and the corresponding voltage Uj were 
displaced about 45° from the line voltage. In order that the voltage supplied to the equalizing circuit's phase- 
sensitive clectronic amplifiers coincide with the line voltage, it was necessary to place a phase-displacing RC 
network at the output of the circuits of Figs. 8 and 9. For taking off the curves of Figs. 8 and 9, one of the 
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potentiometers was turned to the angular position corresponding to a = 0, and the rotor of the second potentio- 
meter was adjusted each time so that a minimum value of the voltage Ures was obtained. As is clear from 

Figs. 8 and 9, the maximum values of Uses were, respectively, 95 millivolts and 145 millivolts. For the parallel- 
connected stator windings, the characteristics of the unbalanced voltage Ug = f(a) in the compensation circuit 
(Fig. 10) were also taken off, This curve is rectilinear throughout the entire working range of the unbalancer 


which, given a linear amplifier, guarantees that the sensitivitiy will be independent of the angular rotation of 
the potentiometer. 
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Fig. 9 


Because of the high accuracy requirements on the linear induction potentiometers, the amplitude and 
phase errors of the prepared potentiometers were measured by means of the compensating circuit shown in Fig. 
11. This circuit is a somewhat simplified version of the compensating circuit of a particular measuring in- 
strument [3]. In the circuit of Fig. 11, the output voltage of the test piece is compensated by the voltage drop 
Ue across the test divider Re, which corresponds to the computed magnitude of the induction potentiometer's 
Output voltage, and by the voltage drop, AU across the linear (wire) divider, Rg. The compensated unbalance 


voltage AU (in amplitude) is taken from divider Rg and, in phase, from the rotor rotation of the phase-shifter 
PS. 


On the basis of the results thus measured, the amplitude error is immediately obtained 


Ry 


AU,% = 100%. 





R 
é max 


The phase error is found from the relationship 


|\A Uelsin 3 
[Ue |\4|8Ujcos* * 





Ag=ranAgo = 


Here, 6 is the angular rotation of the phase-shifter rotor, Remax is the resistance of the test divider cor- 
responding to the compensating current for Uz = Ugmax, Ra is the divider resistance necessary to compensate the 
Output voltage for a definite angular rotation of the induction potentiometer rotor. 


It was demonstrated by the error measurements that, even for relatively crude preparation of the potentio- 
meters, it is possible to guarantee an amplitude error of the first harmonic not exceeding 40.5% of Ugmax in a 
range of 150-160°. The amplitude error decreases with a lowering of the working range of the potentiometer. 


* As in original (cos), Editor. 
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Fig. 10. Dependence of the unbalanced voltage on Fig. 11. Circuit designed to measure potentiometer 
the angle of rotation, for the two parallel-connected error; PS is the phase-shifter, IP is the induction poten- 
potentiometers shown in Fig. 9; Ry = 10 kilohms, tiometer, VM is the vacuum-tube voltmeter, EO is the 
Rg = 140 ohms, Ry = @, electronic oscillograph. 


The phase error (variation of phase of the unloaded potentiometer's output voltage as a function of the 
rotor's angular rotation) did not exceed + 30° even for the worst test pieces. 
The errors cited were affected by technological factors and could be lowered were the potentiometers 
prepared with a higher degree of accuracy. 
SUMMARY 
The following conclusions may be drawn from the results of our investigation, 


1, Adoption of the magnetic system allows for the production of linear induction potentiometers for in- 
dustrial use, satisfying established technological requirements for working angles up to + 90°, with correspond- 
ing linear characteristics. 


2. The computed values for current and voltage for actual potentiometers working in various regimens 
coincide satisfactorily with the experimental data. 


3, The working characteristics of these potentiometers are superior to those of known contact-less position 
transducers currently used in industry with respect to the range in which linearity prevails, symmetry, and smal] 
residual unbalanced voltage in compensating circuits, as affected by the phase error of the potentiometers. 


4. Investigation of the linear induction potentiometers herein described showed their suitability for use 
in automatic control systems, 


5. The use of contact-less induction potentiometers increases the reliability, durability and accuracy of 
industrial automatic control system. 


Appendix 





Example of a Linear Induction Potentiometer Calculation 





Below is given an example of the calculations illustrative of the above method for potentiometers supplied 
by a standard source. 


Given; Uy = 127 volts, Ugnjay = 12 volts, 5Ugmax = + 0.005 and zy = (1000 + j1000) ohms. To be de- 
termined are the basic physical parameters, and the electrical parameters of the potentiometer. 


1. The ratio of the number of turns, according to (5): 





12» / G3 x 100F + 13x 1000 FTX 
c=" 1000 +- 10008 yn 0.008. 
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Here, y is obtained from (4) 


, 0.005 (10002 +4. 10002) 
¥= 26 1000 = 36 





2. Inductive reluctance; 


Y 26 


2° — = 0.12)2 = 1680 ohms 


3. Current in the unloaded potentiometer (7) 
127 


I in “teenie 
1U0L>= 1680 x 2.6 ~~ 0.029 a. 


4. We take Biyax = 2500 gauss, 5 = 0.02 cm, d, /d = 0.8 and d = 5cm., 
5. Magnetizing force of the air gap (8) 
F, =0.8 x 2 x 2500 x 0.8 x 0.02 = 64 a. 


6. Number of turns of the stator winding 


Fy 64 
hy, (0.029 





W,= = 2200. 


7. Magnetic flux linked with the stator winding (9), 


127 x 108 
26x 444 x x2 = 1% mks 





Por = 


8. Block width according to (10)(taking k, = 0.9) 


10* x 2 x 1.25 


4=— mo x0 9-7 cm, @, = 0.35 cm. 


9, Cross-sectional area of the magnetic conductor backs, according to (11) 





10¢ 
So. Seg = 00 = 4 cu. 


10. Area of the copper stator winding, for 4 = 3.5 amp/mm’, 


I 
1UL 0.029 
S1 =~ = Fe = 0.0083 mat, 





We take d, = 0.1 mm, S, = 0.0078 mm’. 


11, Height of one stator according to (12), where we take ke = 0.35, 


0.007810-* x 2200 _ 
hi = 3x09 xX089 TE mM 
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12. Active reactance of the stator coils: 


lw, 4 0,179 x 2200 


r=pP Sy = TF 0.0078 .©~6—Clu 885 ohms 


where 1, » mw (d+ hy) = @ (54 0.7) = 17.9 cm. 
13, Number of turns of the rotor winding 


W, = kw, = 0.125 X 2200 = 274. 


14, The dimensions of the rotor aperture are taken equal to those of the stator aperture 


he = 0.7 cm, c= 0.7 cm. 


15, Area of the rotor winding, for ky = 0.25 


Sinhy =7%5 0.25 
Sy = = aa = 02.0445 ma? 


We take dy = 0.23 mm, S, = 0.0415 mm?, 


16. Active reactance of the rotor coils; 








lw, 1 0.135 x 274 
a=o-F =" FT ta = 16.5 ohms 
where 1, = 9 (d—hg) =m ( 5—0.7) = 13.5 cm 
17, Current in the stator winding (13) for a = a,x 
F 127 ts 
> 27 1680 x 0.1253 ~ 
2(885 +-70.3 x 1680) + 2/1680 [1 — 4033-47 1008 4-2)1680 x 0.425" | 
— 566.5 
where = 0.0272e, wai’ 


5, = 2(16.5 40.3 x 1680 x 0.125%) +1000 + 7 x 1000 = 1033 + 7 x 1008. 


18, Current in the rotor winding (14) for a = a&max 


2/1680 x 0.125 x 0.0272 e~ 366.6 


= — j21° 
ve 216.5 + j0.3 x 1680 x 0.120%) + 1000 + 71000 + 2/1680 x 0.1258 0.00776 * - 





19. Voltage across the load (15) for a= a... 


Oy = 1424) = 0.0077-c— 2" 1415¢%48° — 10.8 pia’, 


_ As is clear from this last expression, the secondary voltage was obtained close to the given one. If neces- 
sary, the magnitude of k may be specified by Equation (5), taking into account the active reactance of the stator 


and rotor windings. 


The computations for the induction potentiometer were carried out by Engineer B. B. Buianovoi. 


272 








LITERATURE CITED 


{1] 0.1. Aven, E. D, Demidenko, S. M. Domanitskii and E. K. Krug, “An electrical servo-mechanism 
with speed regulation,” Automation and Remote Control (USSR) 17, 3 (1956).* 


(2] Iu. M. Pul‘er, “Linear induction potentiometers," Automation and Remote Control (USSR) 17, 7 (1956).* 





[3] lu. M. Pul'er, “Instrument for measuring amplitude and phase errors in ac computing systems," Izmeri- 
tel'‘naia Tekhnika 1956, No. 6. 


Received July 20, 1957 


*See English translation, 


273 





SS 





TELEMETERING COMPENSATION DEVICE FOR LINEAR TRANSFERENCE 


A.A. Kol'tsov and L. F. Kulikovskii 
( Kuibyshev) 


An inductive servosystem for linear transference that may 
be used in automatic and remote control devices is described. A 
brief theoretical analysis and construction data for the system em- 
ployed by the Kuibyshev Industrial Institute are given. 


One of the many possible applications of ferrodynamic measuring mechanisms with independent excita- 
tion and linear displacement of the moving part* is the combined utilization of two such measuring mechanisms, 
forming a telemetering unit for measuring both small and large displacements or other quantities which can be 
transformed into displacements. 
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Fig. 1 





The telemetering unit, represented schematically in Fig. 1, consists of two parts — the master system A 
and the receiving system B. The magnetic circuits of systems 1 and 1° are made of sheet electrical steel. 
Each magnetic circuit consists of an assembly of three rods ~ one central and two side ones (Fig. 2). The ex- 
citation windings 2 and 2' (Fig. 1), fed from the ac supply, are placed in the side rod grooves, Frameless coils 
3 and 3°, which encircle the central rod, are placed in the uniform airgaps of the magnetic circuit. By means 
of special devices these coils are able to move along straight lines. The output terminals of both systems are 
interconnected by line 4, whereas the coils are connected to the terminals by means of light flexible conduc- 
tors 5 and 5° or by sliding contacts. 


The coil of the master system is connected mechanically to the objects whose displacement it is required 
to measure. The coil of the receiving system is supplied with a pointer and the indication on scale 6 gives the 
value of the measured displacement. 


The unit works in the following way. The excitation windings of each system produce magnetic fluxes, 


*The measuring mechanism was proposed and developed by L. F. Kulikovskii and A. A. Kol'tsov, An experi- 


mental model was made at the laboratory of the “Automatic, Remote-Control and Measuring Instruments and 
Devices* Faculty of the Kuibyshev Industrial Institute. 
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equal in amplitude and directed toward each other. Hence magnetic induction in the airgaps of the system 


will remain constant along the whole of their length (Fig. 3), but the induction in the center rod will change 
in the plane of the coils according to a linear law ( Fig. 4): 


B, = kel, (1) 


where kg is the proportionality factor, J is the distance of the plane under consideration from the neutral lines 
0-0 ( 0°-0"). 


The relationships shown in Figs. 3 and 4 are fully 
confirmed experimentally. 








If coil 3 (Fig. 1) is displaced by the measured 
f 7) object through a length 7 from the neutral line 0-0, 
an emf will be induced in the coil 
































A E,; =2ew B,al, (2) 





Fig. 2 where w is the angular supply frequency, wis the number 
of turns in the coil, Bg is the magnetic induction in the 
airgap, a is the dimension of the core (Fig. 2), and! 
is the coordinate of the coil. 


G At the same time an emf will be generated in 
Upper limit coil 3', which remained in the initial position 





P Ey = 2aw'B'a'l’, (3) 





Lower limit as the result of this, a current will flow in the circuit 
. connecting the two coils 








Fig. 3 te (4) 
Z 


where Z is the impedance of this circuit. 


Equality in (4) is achieved by adjusting the excitation current phases of systems 1 and 1° until they coin- 
cide. 


For stable functioning of the system impedance, Z must be inductive, 


Figure 5 is a vector diagram which shows that the interaction between the magneitc induction in the air- 
gaps of the master and receiving systems and component Ig of current I, which is in phase (or antiphase) with 
the magnetic induction, will produce a force acting on the free sides of coils 3 and 3°. 


Force F acting on coil 3° will cause its displacement, which will continue until there arises in the coil an 
emf E} , equal in amplitude but opposite in phase to emf E 1+ The balancing of the two emf's operating in the 
coil circuit leads to the disappearance of current I and the stopping of coil 3’, whose position at J * indicates 
the unknown displacement, since the value of emf E;is proportional to displacement 1‘ ( 3). In other words, 
the coils of the two systems will always assume similar positions with respect to the neutral lines, i.e., | /1'= 
const. This means that the master and receiving systems can in general have different dimensions. The coil 
displacement equation, obtained in the normal way, has the following form: 


dl’ dl’ wr [2Bia'w’|* 


l : 2B, a’w’ 
moe + Pat pe + Fe = 2eoByal 9 


Mis eS (5) 
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where mis the mass of coil 3° and the parts attached to it, p is the damping coef- 
ficient, x is the inductive reactance of the coil circuits, Er is the friction force 


which impedes the movement of coil 3°. 

















5 When coil 3° movement, which can be oscillatory or aperiodic depending on 
. the degree of damping: 

: B= P | 
4B, a’ w’ ( 6) / 
ios se MOI 

Fig. 4 
comes to an end, the coordinates of the receiving coil will be given by the ex- 
pression 
2 , | 
OY me OO 4 —___ fy, | 
B, a'w wz (2B, a’w’ | (7) 


The second term of Equation (7) represents the absolute error of measurement which will only equal zero 
when Ee 0. The last condition cannot be achieved; and it only remains to strive for a minimum friction force, 
whose value is determined by the weight of the coil and attached details, the material and quality of finish of 
the bearing surfaces and those of the guiding wires or sliding contacts. 


The relative nominal error due to friction is given 
by the expression | 
§ 252°F, 
0 ’ 
| — 


t wl y Bi, a’ w’ Boaw 





(8) 


Formula (8) indicates the means of lowering y Fr: 


& E Since the error due to friction is the governing 
—= _ factor in determining the accuracy of measurement, the 
design of the telemetering unit is based on the assump- 
tion that the following quantities are known; the length 

[-—* of the receiving unit scale 2 I'y,, the nominal error of 
measurement K, the coefficient n which is greater than 
one and accounts for errors other than those due to fric- 

5, tion, the friction force Fj, the impedance of the coil 
circuits Z (arg Z = 4/4) and the frequency f of the ex- 
citation voltage. 








Fig. 5 


Formula (9) 


nfl ZF 
Ev = 29.8 V ms (9) 


gives the emf induced in the receiving coil when it is fully deflected. This emf must be equal to the one in- 
duced in the coil of the master system when it is also fully deflected. 


Since the nominal measured displacement / yn is given, the quantity characterizing the master system is 
known to us 


B,aw = on. (10) 
4nfl in 
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as well as that for the receiving system 


(11) 


When the dimensions of the magnetic circuit are determined, the number of turns in the coils and the 


number of ampere turns in the excitation windings do not present any difficulties if the following relationships 
are also used 


c 


B, (8 + A) 


c=a, Ipwp=k, 
F Ho (12) 


where B is maximum permissible value of magnetic induction in the material of the magnetic circuit, o is the 
coefficient of the magnetic flux utilization [1], c is the dimension of the magnetic circuit (Fig. 2), IpWg is the 
number of excitation ampere turns per winding, 5 is the air gap, A is the gap in the butt joints between the 
central and side rods of the magnetic circuit, and k is the coefficient accounting for the reluctance of steel. 


Experiments conducted with a working model of 
this unit demonstrated the high efficiency of the mag- 
netic circuit shown in Fig. 2; the magnetic flux utiliza- 
tion coefficient was equal to 0.7. 


















































Lf? 
hh ' The ac resistance of the coils and their lead-in 
S wires can be considered as purely resistive; hence the 
b. r impedance of the transmission line will be given by 
i5+ \3 the formula 
X : . 
= “ 2 
= on *s% V [fa- e+ r] +3 am 
‘ 4s \ 
10 1 
{ ee INS where ry, and r,, are the resistances of the receiving 
"A “4 \ and master coils and their lead-in wires, respectively. 
x Die Figure 6 shows that there is a good reason for de- 
a5} +z" parting from condition arg Z = y= 2 /4whenr = const., 
* and for increasing the angle y somewhat, since a small 
\ decrease in the controlling force F causes a large de- 
rae a) gdegrees crease in the power P and current I of the coil circuit. 
Fig. 6 This increase in the angle is possible, of course, if the 


resulting impedance modulus does not exceed that of 
Z, which was used for calculations. 


The working model of the telemetric unit consisted of two similar systems (Fig. 7), each of which had the 
following characteristics; 


height of the magnetic circuit 60 mm; 

length of the magnetic circuit 180 mm; 

cross section of the central rod a xX c = 20 x 20 mm?; 
airgaps 6 = 3 mm; 

length of scale 27", = 120 mm; 


magnetic circuit material; sheet electric steel E42; 


excitation winding (each bobbin) had a number of turns Wp = 411; wire used; enamelled, 0.47 
mm diameter; excitation frequency f = 50 cps. 
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induction in the airgap By = 0.08 v sec /m® at excitation current of 0.75 amp; 








moving coil (frameless) had a number of turns W = 500, of enamelled 0.11 mm diameter wire. 


friction force Fy = 0.08 g. 


@ 





0 9 8 


Fig. 7. Telemetering unit receiving system: 1) magnetic circuit, 2) 
excitation bobbin, 3) moving coil, 4) clamping bolt, 5) scale fixing 
bracket, 6) scale, 7) mounting plate, 8) pointer, 9) mounting, 10) 
roller, 11) guiding wire, 12) guide wire fixing clip, 13) guide wire 
tensioning device, 14) core fixing screw, 15) case. 


With a coil-circuit impedance of 270 ohms, the error due to friction was less than 0.5%. The required 
value of angle g was obtained by inserting an inductance into the transmission line. 


LITERATURE CITED 


{1] L. F. Kulikovskii, Electrical Measuring Instruments for Observing Boring Processes (Gostoptekhizdat, 


* Gostoptekhizdat — State Fuel Tech. Press — Publisher's note, 





Received May 31, 1957 








